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Abstract 

The emergent gravity paradigm interprets gravitational field equations as describing the ther¬ 
modynamic limit of the underlying statistical mechanics of microscopic degrees of freedom of the 
spacetime. The connection is established by attributing a heat density Ts to the null surfaces 
where T is the appropriate Davies-Unruh temperature and s is the entropy density. The field 
equations can be obtained from a thermodynamic variational principle which extremises the total 
heat density of all null surfaces. The explicit form of s determines the nature of the theory. We 
explore the consequences of this paradigm for an arbitrary null surface and highlight the thermody¬ 
namic significance of various geometrical quantities. In particular, we show that: (a) A conserved 
current, associated with the time development vector in a natural fashion, has direct thermody¬ 
namic interpretation in all Lanczos-Lovelock models of gravity, (b) One can generalize the notion 
of gravitational momentum, introduced in arXiv 1506.03814 to all Lanczos-Lovelock models of 
gravity such that the conservation of the total momentum leads to the relevant field equations, (c) 
The thermodynamic variational principle which leads to the field equations of gravity can also be 
expressed in terms of the gravitational momentum in all Lanczos-Lovelock models, (d) Three dif¬ 
ferent projections of gravitational momentum related to an arbitrary null surface in the spacetime 
lead to three different equations, all of which have thermodynamic interpretation. The first one 
reduces to a Navier-Stokes equation for the transverse drift velocity. The second can be written 
as a thermodynamic identity TdS = dE + PdV. The third describes the time evolution of the 
null surface in terms of suitably defined surface and bulk degrees of freedom. The implications are 
discussed. 
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1 Introduction 


The dynamical evolution of a fluid or a gas can be studied without any direct reference to the fact that 
they are made of microscopic degrees of freedom, viz., atoms and molecules. Such a description was 
known for centuries before physicists realized that matter is made of discrete entities. But the existence 
of the microscopic degrees of freedom had always left a clear signature even at the macroscopic scales 
in the form of the heat content of matter. While even the cave men knew the distinction between a 
hot body and a cold one, the real nature of heat was not well understood until Boltzmann pointed 
out that the heat content of matter is a direct evidence for the existence of the microscopic degrees of 
freedom. Boltzmann essentially said: “If you can heat it, it must have microstructure”. In other words, 
the microscopic degrees of freedom make their presence felt even at the macroscopic scales (in the form 
of heat) and the correct (thermodynamic) description had taken this into account phenomenologically 
in terms of temperature etc., even before Boltzmann explained to us what it really is. The key new 
element in thermodynamics which is absent in, say, the Newtonian mechanics of point particles, is the 
heat content TS of the matter, which is the difference (F — E) between the free energy and the internal 
energy of the system. In terms of densities, for systems with zero chemical potential which we will be 
interested in, the heat density is Ts = P + p where s is the entropy density, p is the energy density and 
P is pressure. 

Over the decades, it has been realized that spacetimes — through the existence of the null surfaces 
which act as horizons to certain class of observers — also possess the heat density Ts [1-8]. This 
connection between thermodynamics and spacetime dynamics forms the core of the emergent gravity 
paradigm [9-16]. This paradigm rests on the results obtained over the last decade or so which suggest 
that the field equations of gravity, in a large class of theories, have the same status as the equations 
describing, say, a fluid or an elastic solid. That is, gravity emerges in the thermodynamic limit of the 
statistical mechanics of the atoms of spacetime [17-20]. The heat density of spacetime is an evidence for 
the existence of the atoms of spacetime just as the heat density of matter was successfully interpreted 
by Boltzmann as evidence for the atomic structure of matter. There is considerable amount of evidence 
which suggests that this is indeed a useful and correct point of view to pursue [21,22]. 

The development of emergent gravity paradigm has helped us to understand several interesting 
features of classical gravity itself [23], further bolstering our confidence in the veracity of this approach. 
In particular, the following results are of specific interest to the current work. 

1. The first one has to do with the curious relationship between Einstein’s field equations and the 
structure of null surfaces in the spacetime. Previous works have shown that they manifest in 
three different ways: 

• In the most general situation, there arises an identification between Navier-Stokes equation 
and Einstein’s equation. Einstein’s equations, when projected on an arbitrary null surface, in 
any spacetime, leads to Navier-Stokes equation of fluid dynamics [24,25]. (This generalizes 
the previously known results for black hole horizons [26,27].) 

• From the Einstein equations applied to a null surface, one can get [15, 28-31] a thermo¬ 
dynamic identity of the form T6\S = S\E + PS\V in which the symbols have their usual 
meanings and the variation can be interpreted as changes due to virtual displacement of 
the null surface along null geodesics parametrized by the affine parameter A off the sur¬ 
face. Initially proved for a few configurations with a high level of symmetry [31, 32] this 
result has now been generalized for arbitrary null surfaces in both general relativity and 
Lanczos-Lovelock theories of gravity [33,34]. 
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• A comparatively more approximate relationship between the null surfaces and Einstein’s 
equation emerged from the early work [35], which ‘derived’ Einstein’s field equations using 
the local Rindler horizon as a null surface and the Clausius relation. This relies heavily on 
the structure of Raychaudhuri equation as well as the assumptions: (a) the entropy density 
is one quarter of the transverse area and, more importantly, (b) the quadratic terms in the 
Raychaudhuri equation (involving the squares of shear and expansion) can be set to zero. 
Since neither the Raychaudhuri equation nor the assumption that entropy is proportional to 
the horizon area hold for theories more general than Einstein’s gravity, this approach could 
not be generalized in a simple manner to more general class of theories. 

2. The second result pertains to the derivation of gravitational field equations from a thermody¬ 
namic variational principle [14,16]. It turns out that maximizing the sum of gravitational heat 
density and matter heat density, associated with every null surface in the spacetime, leads to the 
appropriate gravitational field equations. This derivation of the field equations is consistent with 
the emergent gravity paradigm and is much more general than the other approaches (like, for 
example, [35]) used in the literature to obtain the same. Not only does the approach has a high 
level of logical simplicity, it also generalizes smoothly to all Lanczos-Lovelock models. 

In this paper, we shall revisit these issues and elaborate further on them using another recent develop¬ 
ment, which was also motivated by the emergent gravity paradigm. This was the introduction of the 
notion of the momentum attributed to the spacetime by a class of observers, proposed in [19,36]. This 
proposal was in the context of Einstein’s theory. But, since virtually every result of emergent gravity 
paradigm could be generalized in a meaningful way to Lanczos-Lovelock models, we would expect the 
notion of gravitational momentum to possess a similar generalization. We will show that this is indeed 
the case. 

Further, it seems reasonable to expect that: (a) The appropriate projections of the gravitational 
momentum must allow us to obtain the above results — describing Einstein’s equations in a thermody¬ 
namic language — in a straightforward manner, (b) The thermodynamic variational principle should 
have a simple representation in terms of the gravitational momentum. 

We will see that these expectations are also borne out. Our analysis will also highlight the thermo¬ 
dynamic significance of several variables which were originally considered purely geometrical. 

To achieve these goals, we will use two different foliations of the spacetime. The first one is the 
standard (1+3) foliation in terms of space-like hypersurfaces determined by constant values of a suitable 
time function t(x). The second is a coordinate system which is adapted to a fiducial null surface. This 
adaptation is closely related to the notion of a local Rindler horizon [35] and will prove to be quite 
useful in illustrating the relationship between field equations and the null surface thermodynamics. We 
will see that the Noether current, gravitational momentum and related constructs, associated with the 
time evolution vector field in the two foliations allow us to obtain the results mentioned above. 

The paper is organized as follows: In Section 2 we introduce the Noether current, gravitational 
momentum and a closely related entity which we call the reduced gravitational momentum in gen¬ 
eral relativity and Lanczos-Lovelock gravity. The next section deals with the two coordinate systems, 
one associated with the (1+3) foliation of the spacetime and the other being the Gaussian null co¬ 
ordinates adapted to an arbitrary null surface. In Section 4 we discuss the Noether current and the 
associated thermodynamic results for both the (1+3) foliation and the Gaussian null coordinates. The 
thermodynamic interpretation of the reduced gravitational momentum is given in Section 5 and the 
thermodynamic variational principle for null surfaces is presented in Section 6. The final section, i.e., 
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Section 7 discusses the three different projections of the gravitational momentum vis-a-vis a null surface 
and their thermodynamic interpretation. We end the paper with a short discussion on our results. 

We will set c, H and 167rG to be unity for most part of our discussion. (Sometimes, when we switch 
to G = 1 units, it will be mentioned specifically.) Thus Einstein’s field equations in this notation takes 
the form 2G ab = T ab . The Latin letters, a,b,c,... runs over all the spacetime indices, the Greek letters, 
/x, i/, a,... runs over all the spatial indices and the upper case Latin letters, A, B,C,... runs over the 
co-dimension two surface. 

2 Noether Current and the Gravitational Momentum 

Given an arbitrary vector field v a in the spacetime, one can define three other geometrical objects 
which depend on it and have direct thermodynamic significance under certain circumstances. The first 
is a conserved current, J a [v\ which can be introduced without invoking any symmetry or invariance 
principle. (This happens to be the usual Noether current; but the usual way of deriving it using 
diffeomorphism invariance of Hilbert action is misleading since it suggests, J“[u] has something to do 
with action and gravitational dynamics. It has nothing to do with either and its conservation is a 
trivial algebraic identity.) The second is the gravitational momentum associated with a vector field 
P a [v] which was introduced in Ref. [36]. The third is a closely related vector (that appears in the 
definitions of both J a and P a ) which we will call the reduced gravitational momentum V a . Each of 
these vector fields, associated with a given vector field v a , can be defined for both Einstein’s gravity as 
well as for Lanczos-Lovelock models. We will now introduce these vector fields. 

2.1 Noether current 

Given an arbitrary vector field v a 1 we can immediately construct a conserved current J a = X7 b J ab from 
the antisymmetric 2-tensor J ab = X/ a v b — X7 b v a . (The normalization of this current is arbitrary; we have 
chosen it so as to make the later results transparent and simple. To match the conventional description, 
the left hand side should be multiplied by 167 tG which we have set to unity.) This construction does 
not require any mention of diffeomorphism invariance or action principles and is completely devoid of 
any dynamical content at this stage. Elementary algebra now leads to the expression: 

V=gJ a (v) = y^[v b (W - vV)] = 2 V=gR a b v b + f bc £vK c (1) 

where we have defined 

f ab = V=gg ab ; = -rs 6 + ± (s c a r d db + s c b r d ad ) (2) 

These variables contain the same amount of information as the metric and the connection but has more 
direct thermodynamic interpretation; see ref. [21]. These expressions are generally covariant because 
the Lie derivative of the connection £ v T^ b , given by 

£ V T a bc = V b V c v a + R a cmb v m (3) 

is generally covariant, making £ v N'' b a generally covariant object. 

There is a natural generalization of this current which can be introduced as follows. We begin by 
noting that the J ab and J a can be expressed in an equivalent form as 

Jab = 2P$VcV d -, J a = 2P$V b V c v d -, P$ = (l/2W a 6t-6i5>) (4) 
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where P“ d b is a tensor (which we will call the entropy tensor for reasons which will become clear later 
on) called determinant tensor. This tensor gives us the Ricci scalar from the curvature tensor: 

R=\ {5 a A - S a d 6 C b )R C a d b = Red Rab (5) 

which shows if we take R = F{R < ^ L ,g i k) to be a function of the (2,2) Riemann tensor R^ and the 
metric tensor < 7 ^, then it is actually independent of the metric tensor. That is, treating and g a b 
to be algebraically independent we can also define P^ d tensor through the following relations: 


pab _ 
*cd 


= (™-\ 

\dR c a d b J 


dR 


dg ab J 


.1 


R cd 


VaP*! = 0 . 


( 6 ) 


Clearly, P£ d has the symmetries of the curvature tensor and is divergence-free in all the indices. 

An obvious generalization of this tensor can be obtained by replacing R by some other arbitrary 
function F(R^,S^) constructed out of and Kronecker delta function, and then define P£ d in an 
analogous manner. That is, we define P® d by the relations 


TDab 
^cd 


(.°L) 

\9R C a d bJ 


Qik 


dF 

dg ab 


= 0 : 


VaP^ = 0 . 


(7) 


This requires finding the most general scalar function F(R^,S^) for which the condition V n P, r !' d = 0 
in Eq. (7) is identically satisfied. This problem can be completely solved [37-39]. It turns out that the 
most general function which satisfies this criterion can be expressed as the sum: 


F — ^ ^ CmFm 

m 


( 8 ) 


where c m s are constants and is given by 


F n 


ra 6 a, 2^2 • • •Q'mbm T>cd t~» c 2^2 TDCrndm 
2 m 0 cdc2d2 .. .Cmdm **'Q'b Jr *'a2b2 ’ ’ ’ ^a m b m 


(9) 


where ^cdc 2 dl cmdm completely antisymmetric m-dimensional determinant tensor. The scalar F 

is constructed out of and the Kronecker delta function 6%, without any metric g ab present in it. 

Then P“ d b for the m-th term in the sum in Eq. ( 8 ) can be obtained directly using Eq. (7) and Eq. (9), 


which leads to, 

pab _ 9F m m ^aba 2 b 2 ...q m fc m pc 2 d 2 p e m d m _ ,^ab 

QR C< ^ 2 m c dC2d2...Cmdm~ r *'CL2b2 ' ’ * "tlm&m ^cd 


( 10 ) 


The mth order term F m can be expressed as F m = Q a cd Rab- f° r m ~ 1j Red an d Qcd coincides. 

Using the definition of p abcd from Eq. (10), in Eq. (4) leads to a natural generalization of J ab and 
J a which, as we shall see later, will be closely related to the Lanczos-Lovelock models of gravity. In 
particular the Noether current in Eq. (4) now becomes 


J a (v) = 2 P^VbVV = 2TZ a h v h + 2P p qra £ v TP r 


( 11 ) 


where TZ% is defined as: Pg = P^R^. 
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The relationship between J ab and J a has an obvious electromagnetic analogy which is useful in some 
calculations to get an intuitive grasp of the results. Given any vector field v a , we can always construct 
the Noether potential J ab (v ) in general relativity, which is exactly analogous to the construction of the 
electromagnetic field tensor F ab starting from the vector potential A> = id (see Eq. (4)). The Noether 
current J a = S7 b .J ab bears the same relation to J ab as the electromagnetic current does to F ab . The 
dual of the tensor J ab is J ab = e abcd J c d and the corresponding (magnetic) current 

J° = S7 b J ab = 2e abcd V b V c v d = e abcd (V fc V c u d - S7 c S7 b v d ) = (12) 

should vanish identically, because this is just electromagnetism in disguise. This is indeed true since 
we have the identity, 

abed f t 1 abed r> , ^ adbc f t , ^ ^aedb rt 

€ J^ibcd — ^ £ -ti'ibcd i ^ £ J^idbc T ^6 -tiicdb 

= —e abcd ( Ribcd + Ridbc + Ricdb ) = 0 (13) 

Further, as in the case of electromagnetism, we can also define an “electric” field and a “magnetic” 
field as measured by an observer with four-velocity u a as 

E a (v\u) = u b J ab (v) = u b ( VV - VV) (14) 

B a (v\u) = \e abcd u b J cd {v) = e abcd u b X7 c v d (15) 

where E a (y\u) stands for “electric” part of J ab [v\ as measured by an observer with four-velocity u a - 
Just as in electromagnetism, we have u a E a (v\u) = 0 = u a B a (v\u) so that E a ,B a are purely spatial. 
The generalization of the electric field to Lanczos-Lovelock gravity is straightforward with p abcd of 
general relativity being replaced by that of Lanczos-Lovelock model. This leads to the expression for 
electric field given by: 

E a (v\u) = u b J ab {v) = 2 P abcd u b V c v d (16) 

The situation with magnetic field is somewhat more complicated in d > 4 dimensions. The J ab has 
d(d — l)/2 components of which the electric vector E a = u b J ab , (with the constraint u a E a = 0) 
contains information about d — 1 components. The magnetic field contains the information about the 
remaining (d — 1 ){d — 2)/2 independent components which is greater than (d — 1) for d > 4. Since 
in Lanczos-Lovelock gravity we are working in d dimensions (d > 4) there is no vector field that can 
account for all the independent components of magnetic field. In d-dimensions the magnetic field is 
given by completely antisymmetric (d — 3) rank tensor B ni -" nd - 3 (v\u) = e ni '" nd - 3abc u a J b c(v). Thus 
total number of independent components is d ~ 1 C d -3 = (d — \){d — 2)/2 as to be expected. Because 
of this rather complicated structure we shall not discuss the magnetic field and shall concentrate on 
electric field, which, fortunately, turns out to be adequate. 

2.2 Gravitational momentum 

The second structure we want to associate with an arbitrary vector field v a is the gravitational four- 
momentum density P a (v ), defined - - in the context of general relativity — as 

P» = -R v a - gV£ v N. £ (17) 
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From Eq. (1) we can substitute for the Lie variation term in Eq. (17), which relates the gravitational 
momentum and the Noether current as 

J a (v) = -P a (v) + 2G a b v b (18) 

The physical meaning of the gravitational momentum can be understood from the following result. 
(This was motivated and discussed in detail in [19,36]. We will not repeat the motivation and logic 
behind this definition here). Consider the special case in which v a is the velocity of an arbitrary 
observer, who will attribute to the matter, with energy momentum tensor T ab , the momentum density 
M a = —T b v b . We would expect the total momentum associated with matter plus gravitation to be 
conserved [36] in nature, for all observers. This condition requires: 

0 = V a (P a + M a ) = V a (- J a + 2G a b v b - T> b ) 

= Va (2 G a b v b - T b a v b ) = (2 G a b - T b ) X a v b = S ab X a v b (19) 

where S ab = (2 G ab — T ab ) is a symmetric tensor and in the last line we have used Bianchi identity 
and the fact that J a and T b are conserved. The above relation should hold for any normalised time 
like vector field v a , which requires S ab = 0, i.e., G ab = 8nT ab , which are the field equation for gravity. 
(This result should be obvious from the fact that X a v b can be chosen to be arbitrary at any event 
even for normalised timelike vector fields. A more formal proof, suggested by S. Date, goes as follows: 
Choose first v a to be a normalized geodesic velocity field with v a v a = — 1 and v a X7 a v b = 0. Then the 
most general S ab which satisfies S ab X a v b = 0 has the form S ab = a(X a v b + X b v a ) + (3v a v b with two 
arbitrary functions a and /3 and an arbitrary vector X a which can be chosen without loss of generality 
to be purely spatial, i.e, v a X a = 0. Choose next the velocity field to be u a = — NX7 a t. Using the form 
of S ab V a u b = 0 leads to a = (3 = 0. This immediately gives S ab = 0.) 

We will now generalise the notion of the gravitational momentum in exact analogy with the way 
we generalised the Noether current, by the relation: 

P a {v) = -TZv a - 2P p qra £ v T p qr 

where 

mn = p;X h = s a b n b a 

Simple algebra gives the equivalent form: 

P a (v) = - J a (v) + 2E£v b 

where J a (v) is the Noether current defined in Eq. (4) and 

K = P%K\ - v/b n ', rnTZ = P c ab R c a d b (23) 

It is possible to prove that E ab is symmetric [40] and that X a E b = 0. To demonstrate the latter, note 
that the covariant derivative of E ab involves two parts, 

V a (. P aijk R bijk ) = P aijk X a R jkbi = P aklm (—V k Ri mab - X b R lmka ) 

= -Vfc ( P kalm R b alm ) + P™XbR r p S q (24) 

d b n = P™V b R r p s q (25) 


( 20 ) 

( 21 ) 

( 22 ) 
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which can be obtained by working in a local inertial frame. These two can be combined to yield 
V a ££ = — Va-Efc , thereby leading to V a E% = 0. 

With this definition of P a , the conservation of total momentum of matter plus gravity leads to, 

0 = V a (P a + M a ) = V a (- J a + 2 E£v b - T£v b ) 

= (2 E£ - T b a ) W a v b (26) 

Invoking the same argument as in the case of general relativity and requiring the above relation to hold 
for all timelike vector fields v a lead to 

K = = l -T b \ (27) 

which are indeed the field equations in Lanczos-Lovelock gravity. Thus imposing the condition that total 
momentum, of matter plus gravity is conserved, with the gravitational momentum given by Eq. (20), 
leads to the gravitational field equations in all Lanczos-Lovelock theories of gravity. This also allows 
us to associate the Noether current in Eq. (11) with the Lanczos-Lovelock models. 


2.3 Reduced gravitational momentum 


We notice that the combinations 


V a = -g ij £ v N^, V a = -2P p qra £ v T p qr ( 28 ) 

appear quite naturally in both Noether current and in the gravitational momentum in Einstein’s theory 
and in Lanczos-Lovelock models. We shall call this combination reduced gravitational momentum V a . 
(We will see later that V a is closely related to the rate of production of heat per unit area on null 
surfaces.) 

The algebraic reason for the occurrence of this combination is as follows. It turns out that, in 
the thermodynamic interpretation of gravity, the combination T ab = T ab — ( 1/2 )Tg ab occurs more 
naturally than the energy momentum tensor T ab with the field equations often arising [ 19 , 20 ] in the 
form 2 R ab = T ab rather than as 2 G a b = T ab in Einstein’s gravity. The total momentum of gravity plus 
matter can be expressed, on-shell , in the form: 

(P“ + M a ) = - g ij £ v N - Rv a - T“u b 

= -g lJ £vN“ + i T8iv b - T b a v b 

= -g ij £ v N? J - = ( V a + M a ) ( 29 ) 

where £4 a = —T b v b is the matter momentum associated with T b . This shows that the vector V a = 
—g l i£ v N£j bears the same relation to T b as P a does with T b . Just as T b appears more naturally in 
the emergent gravity paradigm, the V a will also appear repeatedly in our discussions. 

The notion of reduced gravitational momentum can also be generalized to the Lanczos-Lovelock 
models as well. For pure m- th order Lanczos-Lovelock gravity we have the relation 2[m—(D/2)\L = T. 
So we can write, 


(P a +M a ) = 


qr 


O p qra f pp _. .. 

p v qr 2[m-(D/2)] b 


TS b v b - T?v b 


-2 P p qra 


£ v T p qr - T b V = V a + M° 


(30) 
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where the reduced gravitational momentum is naturally defined as 

V a = -2P p qra £ v T p qr (31) 

In the case of general relativity the gravitational momentum V a defined in Eq. (31) goes over to the 
—g % ^£ v N-j term, as it should. 

3 Coordinate systems and the associated vector fields 

The Noether current and gravitational momentum require a vector field v a for their definition and 
we expect them to have simple physical interpretations when we use naturally defined vector fields 
in the spacetime, associated with the foliations we use to describe the geometry. Our next task is to 
introduce these foliations and the associated vector fields. We will be using two different foliations 
and the corresponding coordinate systems in our work. The first one is based on the standard (1+3) 
foliation while the second one is adapted to a particular null surface. Both of these will turn out to be 
useful in understanding the thermodynamic interpretation of the spacetime dynamics. 

3.1 (1+3) Foliation and the Associated Vector Fields 

This is completely straight forward and we will follow the conventions introduced in [19]. Our primary 
interest is in the case of a coordinate system adapted to a null surface (which we will discuss in the 
next Section 3.2) and the purpose of this section is just to recall the key formulas needed later for 
comparison. 

In the given spacetime, we introduce an arbitrary 1 + (d — 1) foliation based on a time function 
t(x a ) 1 with a unit normal u a (x l ) oc V a t. This splits the metric g a b into the lapse (N), shift (N a ) and 
the ( d — l)-metric h a b = g a b + u a Ub . The unit normal to these hypersurfaces is u a = —lW Q t which 
reduces to — NS ° in the natural coordinate system with t as the time coordinate. Observers with four 
velocity u a will be called fundamental observers. (These observers follow the world lines x a = constant, 
i.e., they have the same spatial coordinates, if we choose a gauge with N a = 0. We will, however, keep 
our discussion general and shall keep N a non-zero unless explicitly mentioned.) This foliation also 
introduces another natural vector field: 

£,a = Nu a —> — (iV 2 , 0); = Nu a -»■ (1, -N a ) (32) 

where the components are in the preferred coordinate system. This vector corresponds to the standard 
timelike Killing vector if the spacetime is static. Further, in any spacetime, it has the time component 
which is unity, i.e., = 1. We will call f a the time evolution vector field. 

The fundamental observers will have (in general, nonzero) acceleration vector = u J V jUi = 
jN/N) which is purely spatial (i.e., u l ai = 0) and has the magnitude a = \foiCd. The conditions 
t{x) = constant, N[x) = constant, taken together, define the codimension-2 surface S with an induced 
metric q a b the area element yjqd d ~ 2 x and the binormal e a i > = r e\ a u^ where r a = e(a a /a) is essentially 
the unit vector along the acceleration. (The factor e = ±1 ensures that the normal r a is pointing 
outwards irrespective of the direction of acceleration. We will usually assume e = 1.) Note that ctj and 
S7iN projected on the t = constant surface coincides. 

We will next consider the construction of an appropriate coordinate system associated with an 
arbitrary null surface which we will use extensively in this paper. 
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3.2 Gaussian Null Coordinates and the associated vector fields 


The construction of a coordinate system associated with an arbitrary null surface has already been 
discussed in detail in [41-43], which we will briefly review. This coordinate system will have the 
following properties: (a) All the redundant gauge degrees of freedom are eliminated, leaving only 6 
free functions in the metric tensor, (b) The null surface we are interested in is chosen to be a surface 
determined by r = 0 where r is one of the spatial coordinates. The other r = constant (but non-zero) 
valued surfaces will represent timelike surfaces and the null surface can be obtained as a limit r —> 0. 
(e.g., this is what will happen in Schwarzschild spacetime if we choose (r-2M) as one of the coordinates 
with the event horizon being the chosen null surface). 

This coordinate system is known as Gaussian null coordinates (henceforth referred to as GNC), 
constructed in analogy with the Gaussian normal coordinates. To handle the fact that the normals 
are null, we need to introduce another auxiliary null vector k a and then construct the coordinates by 
moving away from the null surface along the appropriate null geodesics. After such a construction the 
line element adapted to an arbitrary null surface (identified with r = 0) takes the following form in 
GNC: 


ds 2 = —2 radu 2 + 2 dudr — 2r/3Adudx A + qABdx A dx B 


(33) 


This line element contains six independent functions a, /3a and qAB as advocated, all dependent on 
the coordinates (it,r, x A ). We shall restrict to a > 0, which is adequate for our purposes. The metric 
on the two-surface (i.e. u = constant and r = constant) is represented by qAB- The surface r = 0 is 
the fiducial null surface but surfaces with r = non-zero constant are not null. 

We will now introduce the time development vector £“ appropriate for this coordinate system as 
the one with the components £ a = <5 q in the GNC; that is: 


C = (1,0, 0 , 0 ); G = (— 2m,l ,-rfo) 


(34) 


It can be easily shown that £ a will be identical to the timelike Killing vector corresponding to the 
Rindler time coordinate if we rewrite the standard Rindler metric in the GNC form. Therefore, we can 
think of £ a as a natural generalization of the time development vector corresponding to the Rindler-like 
observers in the GNC; of course, it will not be a Killing vector in general. Since £ 2 = —2ra, we see that, 
in the r —> 0 limit, £ a becomes null. Given £ a , we can construct the four-velocity u a for a comoving 
observer by dividing by its norm \/2rot obtaining: 


u 


l 



Ui = 



1 -rpA \ 
\/2 ra ’ y/2 ra ) 


(35) 


The form of u 1 shows that the comoving observers can also be thought of as observers with (r, x A ) = 
constant. This proves to be convenient for probing the properties of the null surface. 

This four-velocity, has the four-acceleration a 1 = it J Vj u'. The magnitude of the acceleration y/aicd, 
multiplied by the redshift factor \j2ra , has a finite result in the null limit, (i.e., r — > 0 limit): 


N ct| r _>o — 



(36) 


When the acceleration a varies slowly in time (i.e., d u a/a 2 <£/ 1) the second term is negligible and 
Na —> a. The redshifted Unruh-Davies [5, 6] temperature associated with the r = 0 surface, as 
measured by (r, x A ) = constant observer is given [19] by Eq. (36). We will call this temperature as the 
“acceleration temperature”. 
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We will next introduce the relevant null vectors associated with the GNC. Given the four-velocity 
u a and four-acceleration Oi we can construct two null vectors £ l and k l as: 

— \/2vCV — 1 

r = 2---(u a + r Q ); k a = —== (u a — r a ) (37) 

2 \J1ra. 

where r,; is the unit vector in the direction of the acceleration, i.e., r.-, = ai/a. These two vectors £ l and 
k l satisfy: £ 2 = 0, k 2 = 0 and £ a k a = —1 and we have the following components on the null surface: 


r—¥ 0 


r —>-0 


( 0 , 1 , 0 , 0 ) 

/ qABa A a B 

1 ’ 4raa 2 


a-A 



(38) 

(39) 


Since we are essentially interested only in the r —> 0 limit, it is more convenient to work with a simpler 
vector field = Vjr everywhere , which reduces to this £i on the null surface and defines the natural 
null normal to the r = 0 surface as a limiting case. Similarly, we can introduce another vector k a in 
place of k a to simplify the computations. Using the non-uniqueness in the definition of k a , we can 
change it to another vector k a such that, 

ka = ka + Ala + B A&a ( 40 ) 


where are basis vectors on the null surface and l a = £ a - From the property l a e\\ = 0 and £ 2 = 0 
we get £ a k a = —1, since £ a k a = —1. The condition k 2 = 0, leads to a condition between A and Ba 
as: 2 A = qcoB c B D . Choosing A = (qABa A a B /4raa 2 ) and Ba = —{cla/ y/2ra a) leads to the simple 
form k a = (—1,0,0, 0). Thus the vector £ a = V a r and the auxiliary k a = — X7 a u have the following 
components in the GNC: 

£ a = (0,1,0,0), £ a = (l,2ra + r 2 t3 2 ,rP A ) (41a) 

fc„ = (-l, 0,0,0), k a = (0,-1,0,0) (41b) 


Along with these two vectors we also have the vector , which is the time development vector introduced 
earlier. Thus, through the GNC, we have introduced three vectors £ a , k a and £ a . The binormal 
associated with the null surface can be obtained in terms of £ a and k a as e a b = £ a kb — lbk a . 

There are a few more geometrical quantities associated with the null vectors which we will introduce 
for ready reference later on. The first one corresponds to the induced metric q a b on the null surface, 
defined as, 


Qab = gab + £akb + £bk a \ 9b = + £ 0 kb + £bk a 


(42) 


Note that both £ a q£ and k a q£ identically vanishes on the null surface, thanks to the relation £ a k a = — 1; 
we can think of q£ as a projector on to the r = 0 surface, which is two-dimensional. Using this projector, 
we can define extrinsic curvature on a null surface: 

0 ab = C9b V m 4 = \qaQb£tQmn (43) 

If A is the parameter along the null generator £ a on the null surface, we the only nonzero components 
of Q a b are (see Eq. (134) of Appendix A), 


Qab 


1 d 
2d\ 


QAB 


(44) 
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The trace of 0 a b, is 0 = q a b® ab and it is useful to define the trace free shear tensor a ab as, 

&ab — 0a6 ~^Qab® (45) 

Then, as described in Ref. [25] we can introduce shear viscosity coefficient, r] = (l/ 167 r) and the bulk 
viscosity coefficient £ = — (1/ 167 t) as well as the dissipation term by: 

V = 8n (2 Wab a ab + C© 2 ) = 0 a& 0 ab - 0 2 (46) 

The importance of <j a b and T> which will occur repeatedly in our discussion — arises from the 
following fact: It turns out that Einstein’s equations, when projected on to any null surface in any 
spacetime, takes the form of a Navier-Stokes equations [24] with <j ab acting as a viscous tensor and 77 , £ 
acting as bulk and shear viscosity coefficients, In that case, the apparent viscous dissipation is given by 
T>. (The conceptual issues related to this ‘dissipation without dissipation’ since there can be no real, 
irreversible, drain of energy are clarified in [24]). In the GNC, we have on the null surface, 

D = ^q ac q bd d u qabd u qcd - (d u In +<q ) 2 (47) 

which vanishes when d u q a b = 0 on the null surface. We will have occasion to comment on these results 
later on. 


4 Noether current for the time evolution vector field and space- 
time thermodynamics 

I 11 Section 3 we introduced two natural foliations of the spacetime, one based on the (1+3) split and the 
other adapted to a fiducial null surface. These coordinate charts also come with certain natural vector 
fields. In the case of the (1+3) split, the time evolution is related to the vector field £ a we introduced 
in Eq. (32). In the case of the foliation based on a null surface, we again have a natural time evolution 
field £“ given by Eq. (34) as well as the null vector £ a which is tangent to the null congruence defining 
the null surface. It would be interesting to study the Noether current and the gravitational momenta 
corresponding to these vector fields which turn out to have direct thermodynamic significance. In this 
section we shall consider the Noether currents; we will take up the properties of gravitational momenta 
in Section 7. 

4.1 Noether current adapted to (1+3) foliation 

We will begin with a brief description of several properties of the Noether current associated with the 
vector field £ a , developing further the ideas introduced in Ref. [19]. 

Let us start with the situation corresponding to general relativity. The form of the Noether potential 
J ab and current J a is most easily analysed in terms of the electromagnetic analogy introduced earlier. 
Using the foliation vectors u a = —NV a t and £“ = Nu a we readily obtain the following expressions for 
“electric” and “magnetic” fields of J ab [£] as, 

E°(£|u) = u b [V° (Nu b ) - V b {Nu a j\ = -2 Na a - B a {£\u) = 0 (48) 
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For J ab [£] the magnetic field vanishes and the electric field is proportional to the acceleration. So J ab 
can be expressed in terms of E a and u b as 167rJ ab (u) = u a E b (y\u) — u b E a {v\u ) whose divergence gives 
(with 167T inserted; we are now using units with G = 1): 

167rJ a (u) = 16 ttV 6 J ab {v) = u a V h E b (v\u) + E b (v\u)V b u a - E a (v\u)S7 b u b - u b S7 b E a (y\u) (49) 

Projecting along u a leads to the Noether charge density on a t = constant surface given by: 

l6nu a J a (£) = -\7 b E b (£\u) - u a u b \7 a E b (£\u) = ~D b E b (£\u) = D a (2Na a ) (50) 

where Di is the covariant derivative operator on the spatial slices. This matches with earlier results [19] 
and is analogous to V-E = p in electromagnetism. Incidentally, there is a similar result for J a (u); one 
can show that 16nu a J a (u) = D a (Na a ). 

The thermodynamic interpretation of the result 167ro a ,/ a (£) = D a (2 Na a ) can be extracted by 
integrating the Noether charge density u a J a (t;) over a 3-dimensional region V of a t = constant surface 
with integration measure d 3 x\fh to obtain the total Noether charge in that volume. We get: 

f d 3 xVhu a J a (£) = f q a r a d 2 x ( 51 ) 

JV JdV 

where r a is the vector normal to dV and q a is the heat flux vector defined as, 

= (£) (x) < 52 > 

where a“ is the unit vector in the direction of the acceleration. The magnitude of the heat flux vector 
q a gives the heat density, Ts, where T = Na/2it is the local Davies-Unruh temperature [5,6] and 
s = y[q/\ is the entropy density per unit coordinate area. Thus, the total Noether charge contained in 
an arbitrary bulk volume V within t = constant surface is equal to total heat flux contained within the 
boundary dV. When the boundary dV is a iV = constant surface then r a = a a and the total heat flux 
q a r a equals the heat content of the boundary, as derived previously [19]. Also note that 

s/qJ ab {ii) = q a u b - q b u a = 2(Ts)e ab (53) 

where the bi-normal is defined as e ab = (l/2)(d a u b —d b u a )- That is, the Noether potential is proportional 
to the bi-normal of the equipotential surfaces. 

The above analysis uses the fact that u a J a (C) is a 3-divergence, so that the spatial volume integral 
of u a J a (0 can be converted to a surface integral. In this case, it is natural to interpret u a J a ((f) as a 
spatial density, viz., charge per unit volume of space. It turns out that similar results can be obtained 
even for the component of </“(£) in the direction of the normal to the equipotential surface along the 
following lines. It can be easily shown that 

a p J p (0 = - ( g ij - aV) V, (2Nauj) = ■ r/ ; V, (2 Na Uj ) (54) 

where the tensor acts as a projection tensor transverse to the unit vector a 1 . However in order 
to define a surface covariant derivative we need a 1 to foliate the spacetime, which in turn implies 
u l ViN = 0. In this case Eq. (54) can be written as a p J p ( £) = — 'Di(2Nau l ), where Vi is the covariant 
derivative operator corresponding to the induced metric g on the N = constant surfaces with normal 
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a*. (When u l X/iN = 0, we have aj = V ? In IV.) To obtain an integral version of this result, let us 
transform from the original (t, x a ) coordinates to a new coordinate system ( t,N,x A ) using N itself 
as a “radial” coordinate. In this coordinate we have a,; oc 5^ and thus u N = 0, thanks to the 
relation u l ai = 0. Thus T>i(2Nau l ) will transform into 'D a (2Nau a ) 1 where a stands for the set of 
coordinates ( t,x A ) on the N = constant surface. Integrating both sides a P </ p (£) = —V a (2Nau a ), over 
the N = constant surface will now lead to the result (with restoration of 1/167T factor): 



where we have used the standard result \/—g± = N\Jq, with q being the determinant of the two- 
dimensional hypersurface. The right hand side can be thought of as the difference in the heat content 
Q(h) — Q{h) between the two surfaces t = t-2 and t = t\ where: 



This looks very similar to the result we obtained in the case of the integral over u r .J l earlier (see Eq. (51) 
with a a = r“ on the N = constant surface), but there is a difference in the interpretation of the left 
hand side. While UiJ 1 can be thought of as the charge density per unit spatial volume, the quantity 
a p JP(0 represents the flux of Noether current through a time-like surface; therefore, a p J p (£) should be 
thought of as a current per unit area per unit time. We will see later that the flux of Noether current 
through null surfaces leads to a very similar result. 

The generalization of these results to Lanczos-Lovelock gravity is straightforward if we rewrite the 
expression in Eq. (52) for the heat flux q a in the form: 



(When Pf d = (1/2)(5£<5j — 5^8 b c ) the q a becomes parallel to o“, thanks to u a a a = 0 and reduces to the 
previous expression.) To get the results for the Lanczos-Lovelock gravity, we only need to replace P“ d b 
for general relativity with the for Lanczos-Lovelock gravity given in Eq. (10). 

In this case, the antisymmetric Noether potential is constructed by using the entropy tensor Pf b 
introduced through Eq. (6) as, 167rJ ab (u) = 2P c °j’V c u d . Then the “electric” component of J ab {£) turns 
out to be, 


E a {£\u) = 2P abcd u b X7 c t; d = 2P abcd u b [NV c u d + u d V c N} 

= 2P abcd u b u d ( Na c - u c u b \7 b N) - 2NP abcd u b ( K cd + u c a d ) 

= 4 NP abcd u b a c u d = —2Nx° (58) 

where 

= —2Pf d u b a c u d (59) 

Note that u a x a = 0, thanks to the antisymmetry of Pf% in the first two indices making y a another 

spatial vector. (However the “magnetic” component does not vanish in this case and is algebraically 

complicated.) The contraction of </“(£) with u a will again lead to the simple relation u a J a {£,) = 
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— D a E a {£\u) : thanks to the Frobenius identity for u a and complete antisymmetric tensor e abcd . (Of 
course, this is again to be expected from the electromagnetic analogy related to V-E = p.) Thus we 
obtain, by using Eq. (58), the result: 


167ra a J“(0 = D a (2 N X a ) (60) 

Integrating over (d — l)-dimensional t = constant hypersurface with d d ~ x x\[h as the volume measure 
we get, 

J dfi-'xVhUaJ*® = J (^VQ pabc ^ b u c apj r a d d ~ 2 x = j q a r a d d ~ 2 x (61) 

where the heat flux vector q a is now defined as 

«“ = (£) (^W) = (£) (^) X" (62) 

Thus, even in Lanczos-Lovelock gravity we can define a suitable heat flux vector which is spatial (though 
is no longer in the direction of acceleration.) If we consider a region bounded by an equipotential surface, 
then r a = a a and the integrand of Eq. (61) will have the combination P abc ^ a a Ubdpu c which can be 
re-expressed in terms of the bi-normal of the surface. This allows us to express the integrand in the 
form Ts = (Na/2n)s and read-off the entropy density as 

a = -lVQP abcd £a b e cd (63) 


This is in fact an alternate way of defining p abcd which justifies calling it the entropy tensor. 

Using the form of the Noether current it is possible relate the evolution of spacetime to the difference 
between suitably defined surface and bulk degrees of freedom. (This was discussed in detail in ref. [19] 
which we will recall here because we will later obtain a similar result for the null surface.) To do this, 
we start with the relation: 


D a (2Na a ) = 8 t r (2NT ab u a u b ) + u a g bc £^ c (64) 

where we have used the field equation R ab = 87 rT a b, in the expression for u a J a {£) and Eq. (50). The 
Lie variation term [19] is given by (see section A.3 of [19]): 

Vhu a g ij £tN% = —h a b£iP ab (65) 


which depends only on the extrinsic curvature through p ab = Vh(Kh ab — K ab ). For a section of a 
spacelike surface V with boundary dV 1 we define the average boundary temperature as: 


T — _ 

avg A 


d 2 


IdV 


Na 

x ^27 = 


IdV 


d 2 XyfqT\ oc 


( 66 ) 


where T\ oc = Na/2TT is the local (Tolman-redshifted) Davies-Unruh temperature [5,6] and A sur is the 
area of dV. We also define the surface and bulk degrees of freedom by: 


N SUI = / d 2 x^/q = 


IdV 


T 2 

L p 


= = K <67) 


16 







The TV sur counts the number of degrees of freedom on the surface area in Planck units; TVbuik is 
the number of degrees of freedom in a bulk volume if an amount of (Komar) energy Ek omar, is in 
equipartition at the temperature T avg . (We do not assume that the equipartition is reached, of course.) 
Using Eq. (65) in Eq. (64) and integrating it over a bulk volume bounded by TV = constant surface can 
be written in the form [19]: 

~J d 3 xVh h ab £ iP ab = g (TV sur - N h uik) (68) 

which shows that the difference between the surface and the bulk degrees of freedom defined in Eq. (68) 
is what drives the time evolution of the spacetime through the Lie variation of the momentum p ab . It 
will turn out that similar result holds for null surfaces as well. 

We conclude this section with a discussion of the Newtonian limit of general relativity treated 
using the Noether current which has some amusing features. The Newtonian limit is obtained by 
setting TV 2 = 1 + 2 <f>, go a = 0 and g a p = 5 a p, where <f> is the Newtonian potential [12]. Then the 
acceleration of the fundamental observers turn out to be a a = d a (j). Since the spatial section of the 
spacetime is flat, the extrinsic curvature identically vanishes and so does the Lie variation term. Also 
2 T ab u a u b = pKomar = p, which immediately leads to (with G inserted, see Eq. (64)): 

V 2 (j) = AnGp (69) 

the correct Newtonian limit. The same can also be obtained using the four velocity u a . The Noether 
charge associated with u a turns out to have the following expression [19] 

D a a a = 16t Tu a J a (u) = 16nf ab u a u b + u a g bc £ u K c (70) 

In the Newtonian limit the following results hold 2 T ab u a u b = p and u a g bc £ u N bc = ~D a a a (which 
follows from the Newtonian limit of the result u a g bc £^N^ c = ND a a a + 2 a a D a N — Nu a g 1 ^£ u Nij and 
the fact that in spacetime with flat spatial section the term u a g bc £^N^ c identically vanishes). This 
immediately leads to Eq. (69). 

We also see that the Noether charge is positive as long as p > 0 in the Newtonian limit. In fact, 
the Noether charge contained inside any equipotential surface is always a positive definite quantity as 
long as r“ and a a point in the same direction (which happens when T ab u a u b >0). To prove this we 
can integrate the Noether charge over a small region on a t = constant hypersurface to obtain, 

f d 3 xVhu a J a {u) = f d 2 Xy/q2Na a r a = f d 2 x ( (71) 

J inconstant J AT,inconstant J N,t— constant V ^ / \ ^ / 

Since p is positive definite in this case the fundamental observers are accelerating outwards and thus 
r a a a = a. The temperature as measured by these fundamental observers is a positive definite quantity 
and so is the entropy density and hence the positivity of Noether charge follows. 

4.2 Noether current adapted to GNC 

We shall next consider the corresponding thermodynamic interpretation of the Noether current when 
we use the time development vector field adapted to the null surface in the GNC. Let us begin with the 
above result, i.e., Noether charge contained in a bulk region equals the heat content of the boundary, 
which was derived for a subregion of a spacelike surface. It turns out that a similar result holds for a 
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null surface as well. Given the fact that the Noether current corresponding to the time development 
vector led to a nice thermodynamic interpretation, we will consider the object £ a J a (0- (I 11 GNC f a 
has the components in Eq. (34) and l a is given by Eq. (41a) but of course our results are covariant.). 
It then turns out that (see Eq. (187) of Appendix B), 

167 r£ a J a (0 = r (0 = 4= £ (2a Vd) (? 2 ) 


where A is the parameter along the null generator £ a , which in GNC is simply u. This equation can be 
integrated over the null surface with integration measure ^Jqd 2 xdX and leads to, 


where 


d 2 xdX y/qi a J a (Q = / d 2 x ( 


/ ct 


■) (XI 


\2ttJ V 4 


A—A2 


A=Ai 


= Q(A 2 ) - Q(Ai) 


« A) = / A (S) ( f)=f ^ T » 


(73) 


(74) 


is the heat content of the null surface at a given A. 

This again shows that total Noether charge density for the vector field f a integrated over the null 
surface equals the difference of the heat content Q of the two dimensional boundaries located at A = A 2 
and A = Ai. Previously the connection between bulk Noether charge to surface heat density was 
derived in the context of spacelike surfaces. The result in Eq. (73) generalizes the previous connection 

between bulk Noether charge and surface heat density in the context of spacelike surfaces - by 
showing that the total Noether charge on a null surface is also expressible as the heat content of the 
boundary. 

The following aspect of this result is worth highlighting. We obtained earlier two results (see Eq. (51) 
and Eq. (55)) of similar nature. The first one (in Eq. (51)) was for a spatial region V contained in a 
space-like hypersurface. In that case, the normal to the surface was u a and the natural integration 
measure, for integrating the normal component of a vector field is u a Vhd 3 x. We computed the integral 
J a u a \fhd?x and found that it is given by a boundary term; we could have also computed the integral 
in a region contained within two boundary surfaces like, for example, in the shell-like region between 
two spherical surfaces of radii R\ and i? 2 . We would have then found that the Noether charge in the 
bulk region is the difference between the heat contents of the two surfaces. 

In the second case, (in Eq. (55)) we considered the flux of Noether current through a timelike surface 
with normal a 1 . In this case we calculated the integral of d p J p {£f) with the measure d 2 xdt^/—g± on a 
timclikc surface and got a similar result. We also mentioned that in this case, the integrand a p J p (Jf) is 
to be thought of as heat flux per unit time. 

In the case of a null surface, our result is similar to the second one, given by Eq. (55). Now 
the corresponding integration measure for integrating a vector held over a null surface is given by 
(, a y/qd 2 xd\ where £ a = dx a /d\ is the tangent vector to the null congruence defining the null surface. 
Therefore, in this case, we calculate the integral over J a i a ^/qd 2 x dX. This leads to the difference A Q 
of the heat content at the two boundaries corresponding to A = Ai and A = A 2 . In the integrand in 
this case, J a l ay /qd 2 x dX one of the coordinates A is similar to a time coordinate rather than a spatial 
coordinate. So we cannot think of J a £ a as charge per unit volume ; instead it represents charge per 
unit area (flux) per unit time and more appropriately, it is the rate of production of heat per unit area 
of the null surface. 
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It is possible to proceed further and relate the change in the heat content with the matter energy 
flux through the null surface. To do this, we use Eq. (1) and field equation G ab = 87 xT ab to obtain (on 
the null surface): 


16nT ab re b = 16 7r£ a J a (0 - l a g bc £iK c 


(75) 


The Lie derivative term can be computed directly to give (see Eq. (225) of Appendix B): 

y/ql a g ij £tN§ = -q ab £z U ab + ^ (76) 

where LP b = y/q[Q ab — g ab (0 + k)] is the momentum conjugate to q ab and 0 ab = q^q^'V m l n where 0 is 
the trace of 0 ab . Integrating this result over the null surface between A = Ai and A = Ai, and assuming 
for simplicity that the boundary terms at A = Ai, Ai do not contribute (which assumes dA/dX = 0 at 
the end points where A is the area of the null surface), we get: 

- Y^ / d 2 xd\ q ab £zU ab = [Q{ A 2 ) - Q(Ai)] - J d\d 2 xJqT ab tt b (77) 


This expression shows that the evolution of the spacetime, which is encoded by the evolution of the 
momentum II ab is driven by the difference between (i) heat content Q at the boundary and (ii) the 
matter heat flux flowing into the null surface. We can rewrite this in a nicer manner as follows. We 
define the surface degrees of freedom as (in units with G = 1 ): 


and the average temperature as: 



We also introduce the effective bulk degrees of freedom by: 


(78) 


(79) 


-Wulk 


1 

(l/2)T avg 


d\d 2 Xy/q2T ab £ a l b 


(80) 


If the matter heat flux, given by the integral in the right hand side, thermalizes at the average temper¬ 
ature of the null surface then, IVbuik will represent the the effective equipartition degrees of freedom. 
We now rewrite Eq. (77) as, 


__L y cf xd \ q ab£i n“ b = 1 y d 2 x (£) VS - i y d 2 x (^) V9 - 2 j d\d 2 x^T ab r( b 

which, on using our definitions, becomes: 


1 

87T 


d 2 xd\ (q ab £ i U ab ) 




-Wmlk 


(81) 


(82) 


This has the interpretation that and hence the evolution of spacetime on a null surface, encoded in the 
Lie variation of the momentum II ab along the time development vector, can be thought of as due to the 
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difference between surface degrees of freedom and the bulk degrees of freedom. This an exact analogy 
to the corresponding result in (1 + 3) foliation presented in Eq. (68), which was originally obtained 
in [19]. 

For the sake of completeness, we clarify the notion of T ab £ a £ b being the heat flux through the null 
surface. This concept has been introduced in [35] and arises as follows: Let there be a matter field, 
in the spacetime, with energy momentum tensor T ab . Around any given spacetime event V, we can 
construct local inertial and hence local Rindler frames. We then have an approximate Killing vector 
field generating boosts, which coincides with the null normal £ a at the null surface (see Section 3.2). 
The heat flow vector is defined as the boost energy current obtained by projecting T ab along £ b yielding 
T a b^ b - Thus the energy (heat) flux through the null surface will be: 

Q = J (T ob £ b ) d,Z a = J T ab i b r^d 2 xd\ (83) 


where yjqd 2 xd\ is the integration measure on a null surface generated by null vectors £ a , parametrized 
by A. Hence, in the null limit, T ab £ a £ b (when £ a —> £ a on the null surface) represents the heat flux 
through the null surface. 

The same argument can also be presented along the following lines. On a null surface we can 
decompose T ab £, b in canonical null basis as, T a b£ b = A£ a + Bk a + CAe-a ■ Then the heat flux through 
the surface is given by the component B along k a , which is off the null surface. This component B is 
obtained by contracting with £ a (since £ 2 = 0 and £ a k a = — 1). This leads to the heat flux through the 
null surface to be T ab £ a £ b . 

5 Reduced Gravitational momentum associated with the time 
development vector 

The expression for the Noether current contains a Lie derivative term which was defined earlier as the 
reduced gravitational momentum. In this section, we shall describe some key properties of this reduced 
gravitational momentum vector in different contexts, emphasising the results in GNC. 

In the case of (1+3) foliation using spacelike surfaces with normal u a = —NV a t (N is the lapse 
function) and induced metric h ab = g ab + u a Ub the reduced gravitational momentum can be related [19] 
to the Lie variation of p ab = \fh{Kh ab — K ab ) by: 

-Vhu a v a {0 = ^ ll „g'-'£,.\;' l = - h ab £ iP ab (84) 

In the case of null surfaces, one can obtain a similar relation. For a general, non-affine parametrization, 
i.e., when the null generator £ a satisfy the relation £ b X/ b £ a = n£ a , we find that the Lie variation term 
is given by (see Eq. (225) of Appendix B): 


-y?4^ a (0 = VqLg ij £^ = -<za 6 ^n ab + ^ (85) 

where n ab = y/q[Q ab — g ab (0+fc)] is the momentum conjugate to q ab - Thus as in the case of the spacelike 
surface, for null surfaces as well, the Lie variation of A+ b is directly related to the Lie variation of the 
momentum conjugate to the induced metric on the null surface. But in the case of null surfaces there 
is an extra term which contributes only at the boundaries A = Ai, A 2 . 
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It can be seen from straightforward algebra (see Eq. (163) of Appendix A) that q a b£^P ab is directly 
related to the object V = 0 a f,0 ab — 0 2 defined in Eq. (46) which, as we mentioned earlier, can be 
interpreted as dissipation [25]. This leads to an alternative expression on the null surface for the Lie 
variation term in the adapted coordinate system as (see Eq. (210) and Eq. (224) in Appendix B), 

Lg ij £^ = 2V+ 4 = dl (y/q) + 2 d x a (86) 

Integrating this expression over the r = 0 null surface with integration measure drxdu y/q, neglecting 
total divergence and dividing by 167 t leads to, 

J d 2 xd\ y/q£ a g lJ £^N" j = j-j d 2 xdXy/qV + J d 2 x d (£) (87) 

which explicitly shows that the Lie variation term integrated over the null surface, leads to the dissipa¬ 
tion and the sdT term (interpreted as dT = (dT/dX)dX). Hence the reduced gravitational momentum 
on the null surface, can be given a natural thermodynamic interpretation. 

What is important regarding the above result is that for d\ q a b = 0 (i.e., the induced metric on the 
null surface is independent of the parameter along the null generator) the dissipation term vanishes 
and thus Eq. (87) can be written as: 

y^ J d 2 xdX y/q£ a g ij £^ = Q(X 2 ) - Q( A,) (88) 

Hence the Lie variation term in this particular situation is equal to difference in the heat content 
allowing us to relate V a to the rate of heating of the null surface. 

6 The thermodynamic variational principle for the field equa¬ 
tions 

It has been shown earlier [14,16] that the gravitational held equations can be obtained by extremising 
the total heat density of all the null surfaces in the spacetime. The purpose of this section and the next 
is to show that this variational principle takes a simple form in terms of the gravitational momentum 
for both general relativity and the Lanczos-Lovelock models. 

6.1 General relativity 

In Section 2 we have defined the gravitational momentum P a (y ) and the matter momentum A4 a (v) 
associated with the vector held v a . Using these two we can construct a thermodynamic variational 
principle leading to the held equations for gravity. We start by dehning a suitable expression for the 
heat density associated with a null surface as: 

Q = J d 2 xdXy/q [{-i a P a (£)} + {-£ a M a {£)}\ (89) 

where £ a is the null vector, such that on the null surface £ a £ a = 0. Then from Eq. (18) we observe that 
the Ricci scalar does not contribute since £ a £ a = 0 on the null surface. From the dehnition of matter 
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momentum as A4 a (£) = —T££ b , we find that: 

Q = J d 2 xdX^/q (t, ab tl h + ^ i a g hc £iN^ 

To understand the interpretation of this quantity as the heat content of the null surface, note that 
the matter term involves an integral over d 2 xdX^/q £ a (T£ £ a ) on the null surface. As we argued earlier, 
T ab £ a £ b can be related to be heat flux due to the matter and hence this integral represents the heat 
content of the null surface contributed by matter. (In the case of an ideal fluid in GNC, T ab £ a £ b will 
be equal to p + p, which is same as the heat density Ts when we use the Gibbs relation.) Similarly, 
the expression obtained by replacing A4 a by the gravitational momentum P a can be thought of as the 
heat density due to gravity. The explicit form involving reduced gravitational momentum supports this 
interpretation. Since £ a = V a </>, we have J a {£) = 0. Hence we have £ a g bc £(.N bc = —2R ab £ a £ b , which 
allows us to write: 


e = /^ 

Varying the integrand with respect to all the null vectors fields £ a , after adding a Lagrange multiplier 
term X(x)£ a £ a to enforce the condition £ 2 = 0 will lead to the field equation G a b = S>TrT ab + Ag ab where 
the cosmological constant A appears as an integration constant. (The details of this derivation has 
been given in several previous works [14,16,19] and hence is not repeated here.) So, starting from the 
projection of gravitational momentum and matter momentum along the null normal i a , we can derive 
the held equations for gravity. 

The same result arises even if we have started using the contractions —£ a P“(£) and —£ a A4 a (£). As 
we have shown in Eq. (149) in Appendix A the Lie variation term and R a b£ a £ b differs only by a total 
divergence. Hence after neglecting the surface contributions, the action would be identical to Eq. (91) 
and thus on variation of £ a it would lead to the held equations for gravity. 

6.2 Lanczos-Lovelock gravity 

The power of this analysis becomes apparent when we realise that the same expression as presented 
in Eq. (89) leads to the Lanczos-Lovelock held equations when we use the corresponding momentum 
P a (£) given by Eq. (20). Using this, we can construct a variational principle associated with a null 
surface, leading to the held equations for Lanczos-Lovelock theories of gravity. This is achieved by 
using exactly the same expression for the heat content: 

Q = / d D ~ 2 xdX ^ \{-£ a P a {£)} + {-£ a M a {£)}} (92) 

where £ a is the null vector, such that on the null surface £ a £ a = 0. As in the case of general relativity, 
the scalar 1Z present in gravitational momentum does not contribute since £ a £ a = 0 on the null surface. 
From the dehnition of matter momentum as A4 a (£) = —T££ b , we arrive at the explicit form: 

q = J d D ~ 2 xdx ^ (r ab £ a £ b + ^£ a p r ^£ e r;^ ( 93 ) 


xdX ( T ab £ a t - ^Rab£ a t 


(91) 
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which gives the variational 


Since £ a = V a </>, we have J a {£) = 0. Hence £ a P r ^ qa £ e T^ q = - 2TZ ab £ a £\ 
principle for the null surface to be: 

Q = J d D ~ 2 

Again varying the integrand with respect to all the null vectors fields £ a , with a Lagrange multiplier 
X£ a £ a to impose the constraint £ 2 = 0, we obtain the field equation to be E ab = 8 ttT ab + A g ab as before. 
Hence starting from the gravitational momentum and matter momentum along the normal, null vector 
field £ a , we can derive the held equations for gravity. 

Thus, one can write down a thermodynamic variational principle directly in terms of the gravita¬ 
tional momentum. The logical simplicity of this result and the fact that it holds for Lanczos-Lovclock 
models without any modification is note worthy. Because of this feature, this procedure seems to be 
the natural way of obtaining the held equations in the emergent gravity paradigm. 

7 Projections of Gravitational Momentum on the Null Surface 

We shall now take up a further application of the concept of gravitational momentum. Given the 
thermodynamic signihcance of the null surfaces, we would expect the how of gravitational momentum 
vis-a-vis a given null surface to be of some importance. To explore this, we have to hrst choose a suitable 
vector held using which the gravitational momentum can be dehned. The most natural choice — as 
before — is the time evolution held £ a . Further, in the canonical null basis (£ a , k a , e^) the gravitational 
momentum can be decomposed as: P a = A£ a + Bk a + C A e\. These components A, B and C A are 
related to the three projections of P a by A = —P a (£ t )k a , B = —P a {£)£ a and C A = P a (t;)e A . So we 
need to consider the following three components q b P b {Q,k a P a {t;) and £ a P a (0 to get the complete 
picture. We will see that each of them lead to interesting thermodynamic interpretation. In view of 
the rather involved calculations, we will hrst provide a summary of the thermodynamic interpretations 
of these projections: 

• The component q b P a (0 leads to the Navier-Stokes equation for huid dynamics, using which 
we can obtain yet another justihcation for the dissipation term introduced in Eq. (46). This is 
described in Section 7.1 

• The component fc a P a (£), when evaluated on an arbitrary null surface leads to a result which can 
be stated in the form: TdS = cLE + PdV , i.e., as a thermodynamic identity. This helps us to 
identify a notion of energy associated with an arbitrary null surface. We obtain this result in 
Section 7.2. 

• Finally, the component £ a P a {Q yields the evolution of null surface, which involves both ds/dX 
and dT/dX , where s is the entropy density and T is the temperature associated with an arbitrary 
null surface and A is the parameter along the null generator £ a . This is studied in Section 7.3. 

We shall now show how all these three results tie up with the notion of gravitational momentum 
and arise naturally from the three different projections of the gravitational momentum. 


xdX^q ( T ab £ a £° - — 


(94) 
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7.1 Navier-Stokes Equation 

The equivalence of field equations for gravity, when projected on a null surface, and the Navier-Stokes 
equation is an important result in the emergent paradigm of gravity. This result, which generalizes the 
previous one [27] in the context of black holes, shows that when Einstein’s equations are projected on 
any null surface and viewed in local inertial frame, they become identical to the Navier-Stokes equation 
of fluid dynamics [24]. Here we will project the gravitational momentum P a on an arbitrary null surface 
parametrized by GNC and show that the Navier-Stokes equation is obtained. For this purpose we will 
consistently use the vector field £ a . Even though we will present the results for the adapted coordinates 
to a null surface (i.e., in GNC), the same result continues to hold in any other parametrization as well. 

In order to project the gravitational momenta 167 tP c (£) = g ab £^N^ b + R£ c on the null surface 
we need to determine the induced metric q bl which [43] turns out to be diag(0,0,1,1). Hence it is 
the angular part which is going to contribute. The vectors i a and k a are already given in Eq. (41a) 
and Eq. (41b) respectively. In the Navier-Stokes equation two other vectors play a crucial role. These 
vectors and their components in the adapted GNC system has the following expressions: 

P (95) 

and 

tta = Wa + aka = ^0, 0, 

From [24] we know that /3a can be interpreted as the transverse velocity of observers on the null surface. 
(In particular it can be interpreted as velocity drift of local Rindler observers parallel to the Rindler 
horizon.) With this physical motivation, let us now start calculating the projection of P a on the null 
surface. Using the coordinates adapted to a given null surface, we arrive at: 

-167T qZP b (Z) = q b g pq £$N b q = 2£ ( N* + q BC £^N A C (97) 

where we have used the identity £ a q b = 0. We next need to find the Lie variation of N£ b . For that we 
use the expression in Eq. (3) for the Lie variation of Christoffel symbol with respect to an arbitrary 
vector field v a and the result: N bc = Q b d T e cd + Qc d £ bd - Then the Lie variation of N bc becomes 

n ]\ra \ad n pe ■ r\ad n pe 

X'vWbc ~ V be c d ~r V ce ot '^ i bd 

= l (<^V c V rf u d + 6 a c \7 b \7 d v d ) - 1 (V 6 V c u“ + V c VhU a ) - i (. R\ mc + R a cmb ) v m (98) 

We obtain the expressions for Lie variation of N^ r . and N^ c on the null surface (located at r = 0) to 
be (see Eq. (191) of Appendix B) 

£eN* = \d u p A = \q AB d u p B + \p B d u q AB 

£^ B g = + 2 S A d B Q — d u T BC 

Substituting these results in Eq. (97) we arrive at, 

-167 Tq p P b = qP (g bc £zK c ) = q PB d u p B + f3 B d u q PB + q PC d C & ~ q BC d u f P c (101) 


(99) 

( 100 ) 


2 13 A 


= ( 0 , 0 , 


(96) 


= £ m S7 m k b = 


a,0, -Pa 


= 0,a, 
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In order to get the Navier-Stokes equation in its familiar form we need to lower the free index in 
Eq. (101) and multiply both sides by (—1/2). Using Noether current for f a we have from Eq. (1) 

qab{g pq £^N b pq ) = q ab J b (f) - 2 R mn r q n a (102) 

On using Einstein’s equations R a b = 8n(T ab — (1/2 )Tg ab ) and the result i a q% = 0, Eq. (101) leads to 
the following result: 

Bt TT mn e m q™ = i/3 A d u In y/q + \d u P A - d A a - ^d A d u In yjq + \ > q A Bq PQ d u tp Q 
= In yfq + -d u p A - d A a - -d A d u In yfq 

- i q A B£p Q d u q PQ + ( q CD d u q AC ) + i d u q CF d c q A F - ^ qAB^u ( q BD d D In yfq) 

= i PaO u In y/q+ i. d u /3 A - d A a - d A d u In yfq- i q AB d u q BD ( d D In y/q) 

+ \ d D ( q CD d u q AC ) + i d u q CF d c q A F - i d c q A Fd u q CF - ^q CD T AC d u qED (103) 

where in the second line we have used the relation Eq. (175b) and in the third line Eq. (175c) as 
presented in Appendix B. From Eq. (174) of Appendix B we obtain: 

^d-upA + d B (^q BC d u q AC ^j + i q CD d u q AD d c In y/q - ^q BD d u q CD £ AB 

+ d u In yfq-[3 a - d A d u In yfq - d A a 

= q™£ e n n +D m a™ + Gn a -D a (j + a'j (104) 

After some trivial manipulations in Eq. (103) and using Eq. (104) the following final expression is 
obtained: 

Btt T mn t m q n a = qf£ t n n + D rn off + 0U a - D a + a) (105) 

which can be interpreted as the Navier-Stokes equation for fluid dynamics. 

The correspondence between Eq. (104) and Eq. (105) with Navier-Stokes equation for fluid dynamics 
is based on the following identifications of various geometric quantities on the null surface: (i) The 
momentum density is given by — fl a /8n. In the coordinates adapted to the null surface, f l a has only 
transverse components which are given by (l/2)p A . This reinforces our interpretation of p A as the 
transverse fluid velocity. Moreover, we can identify (ii) the pressure k/8tt, (iii)the shear tensor defined 
as a mn (see Eq. (45)), (iv) the shear viscosity coefficient r] = (1/167 t), (v) the bulk viscosity coefficient 
f = —1/107T and finally (vi) an external force given by F a = T ma l m . Thus Eq. (104) has the form of a 
Navier-Stokes equation for a fluid with the convective derivative replaced by the Lie derivative. Since 
this equation and its interpretation have been extensively discussed in the references cited earlier, we 
will not repeat them and will confine ourselves to highlighting the dissipation term. 

In order to interpret the dissipation term we start from the heat density q = —4P//’V„£ c V b £ d 
(i.e., heat content per unit null surface volume yfqdfxdX) used in the variational principle, where 
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P$ = (1/2 )(S-S b d - 5 a d 8 b ) and i a is the null generator of the null surface. To connect up the heat 
density with the dissipation term [25] we consider a virtual displacement of the null surface in which 
the volume changes by <5 AAA, where 8 A is a small area element on the two-surface (i.e., qcPx ). Then 
qSAdX represents the change in heat content of the null surface, which is obtained by multiplying the 
heat density q with the infinitesimal 3-volume element on the null surface. Expanding the expression 
for the heat density and introducing temperature through local Rindler horizon [25] we obtain, 

= —SAdX {2 Wab a ab + C0 2 ) + ~dSA (106) 

where the first term represents the (virtual) dissipation of a viscous fluid during the evolution of the 
small area element 8A of the null surface along the null generator and has the following expression, 

dE = 8AdX (2r]<T ab cr ab + C© 2 ) = —SAdX V (107) 

87T 

So the combination V = (0 a {j© ab — 0 2 ) represents the dissipation as we have mentioned earlier. The 
second term in Eq. (106) can be interpreted as (1/2 )kT dN , where T = (k/27t) and dN is the degrees 
of freedom on the null surface corresponding to the change in area 8A. For affinely parametrized null 
generator k = 0, which would lead to, 

= —8AdX (2 V a ab a ab + C© 2 ) = ^-SAdX V (108) 

o7T v ' o7T 

Thus for affinely parametrized null generators the change in the heat content is solely due to the 
dissipation term T>. 

Having derived Navier-Stokes equation from the projection of gravitational momentum P a (£) on 
the null surface, we will now take up the task of projecting it along k a and l a respectively and retrieve 
the thermodynamic information encoded in them. 

7.2 A Thermodynamic Identity for the null surface 

It has been shown that in a wide class of gravity theories, the gravitational field equations near a 
horizon imply a thermodynamic identity TS^S = S\E + P8\V where the variations are interpreted 
as being due to virtual displacement of the null surface along the affine parameter A of k a . This 
result was first obtained for general relativity and Lanczos-Lovelock theories of gravity when (i) the 
spacetime admits some symmetry, e.g. staticity or spherical symmetry and (ii) has a horizon. This 
may suggest that this connection — between the field equations and a thermodynamic identity — is 
a specific phenomenon that occurs only in solutions containing horizons. But this illusion is broken 
in [33] for general relativity and in [34] for Lanczos-Lovelock theories of gravity. There it was shown 
that gravitational field equations near any generic null surface in both general relativity and Lanczos- 
Lovelock theories of gravity lead to a relation: TS^S = 5\E + P8^V. 

Here we will show that this thermodynamic identity is also contained in the gravitational momentum 
for £“ and can be retrieved from its component along £ a which is contained in the projection P a (^)k a 
(which picks out the component along £ a because k 2 = 0 and £ a k a = —1). Since all the details are 
similar to the one in [33] we will be quite brief and just indicate the manner in which the result can be 
obtained. (For detailed derivation see Appendix C of [33]). From Eq. (18) we obtain, 

-k a P a (0 = k a J a (0 - 2GU b k a (109) 
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In the second term we can use the field equations, i.e., 2G a b = T a b leading to the combination T a b^ a k b , 
which is the work function (effective pressure) for the matter (see e.g., [44,45]) which, when integrated 
over the two-surface, will yield the average force F due to matter flux on the surface. The first term, 
viz., the projection of the Noether current k a J a {£) has been evaluated explicitly in the earlier work [33]. 
Using this result from ref. [33] for the projection of the Noether current we arrive at: 

F5X = TS x S - S- X E , (110) 

In this expression, F stands for the average force on the null surface, T corresponds to the null surface 
temperature obtained using local Rindler observers, S = (A/ 4) is the entropy and E represents the 
energy given by: 

E =\lHi)~kl i J dX I■ an) 

where x stands for the Euler characteristics of the two-surface. (For a detailed discussion see [33]). 
A simpler expression covering most of the interesting cases is obtained by setting (a) /Ut|r=o = 0 and 
(ii) 8aol = 0 on the null surface [33]. On imposing these conditions we arrive at the following simpler 
expression for energy as, 

E = \! dX (-2)-kI=\S iX (I) - k 9 ' A < 112 > 

Hence the projection of the gravitational momentum along i a is equivalent to thermodynamic identity. 

7.3 Evolution of the null surface 

Finally, we consider the component of the gravitational momentum P a [£] along k a . If we expand any 
vector in the basis (£ a , k a , e A ) as v a = Ai a + Bk a + C A e A , the component along k a , is obtained by the 
contraction £ a v a , since £ 2 = 0 and l a k a = —1 on the null surface. For a given null surface we will use 
our adapted coordinates, i.e., GNC and will show that this component is intimately connected with 
spacetime evolution. 

In the adapted coordinate system the scalar — f Q P a (£) has the following expression, 

-167T4P a (£) = Lg bc £iK c = 2 £^N r ur + q AB £ i N r AB (113) 

where the first line follows from the fact that on the null surface £ 2 = 0. Both the Lie variation terms 
have been calculated in Appendix B explicitly. From the expressions obtained there the projection of 
gravitational momentum turns out to be (see Eq. (210) and Eq. (224) of Appendix B), 

-16t r4P“(0 = Lg bc £^ c = 2d x a + 2V + -L d 2 xx /q (114) 

v? 

Here T> represents the dissipation term obtained through the Navier-Stokes equation and has the 
definition V = (0 Q b0 ab — 0 2 )- 

To get a physical interpretation we integrate this expression over the null surface with volume 
measure y/qd 2 xd\ (for the null vector £ a the parameter A is just u) and divide by proper factors of 7r, 
and ignore the surface term (which arises from the third term) to obtain: 

J dXd 2 x^l a g ij £ t :N°- j = J dXd 2 x^V + j d 2 x sdT (115) 
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where s = y/q/A and cLT = (dT / dX)dX. This leads to the result that l a P a (0 and £ a g bc £^N^ c can be 
interpreted as the heating rate of the null surface per unit area. Integrated over the affine parameter 
dX and the proper area y/qd 2 x, it leads to the heating due to the dissipation (given by the first term 
in the right hand side) and the integral of sdT which is the second term. 

Incidentally, a similar interpretation can be given for the matter energy flux which crosses the null 
surface as well. Using a corresponding expression for R ab £ a £, b and using the field equations, we can 
easily show that (see Eq. (211) of Appendix B) 


T ab re b 



1 11 d 2 Jq 

— V - — 

87 t 87 t yfq du 2 


(116) 


Integrating both sides over the null surface and ignoring boundary contributions at the ends of inte¬ 
gration of affine parameter, we get 


J d 2 xdX y/qT ab £ a £ b + J d 2 xdXy/qV = J d 2 xTds 


(117) 


This tells us that the heating due to matter flux plus the heat generated by the dissipation is equal to 
the integral of Td u s over the null surface. This reconfirms the earlier interpretation of the projection 
of the momentum contributing to the heating of the null surface. 

These results can be used to re-express the heat content of the null surface which was used in the 
thermodynamic variational principle. Two equivalent forms of the variational principle, which differ 
by a total divergence can be given based on R ab £ a l b and the Lie variation term. These two variational 
principles (neglecting surface contributions) have the following expressions 


Qi = j dXdrxyfq ( -—R ab £ a t + T ab t£ b ) = / dXd 2 Xyfq 


^-V + T ab £ a £ b 

L87T 


and, 

<?2 


JdXd 2 Xy/q 


t a g i3 £tN% + T ab £ a £ b 


16n 


= / dXd 2 Xy/q 


-8ir 


v + T ab n b 


- / d 2 x Tds (118) 


+ / d 2 x sdT (119) 


Note that both these variational principles have the dissipation term T> and matter energy flux through 
the null surface T ab £ a £ b in common. However Q i is connected to TdS while Q 2 is connected to SdT. 
Thus both the variational principles have thermodynamic interpretation. 


8 Conclusions 

Since we have described the physical consequences of the results in the various sections themselves, we 
shall limit ourselves to summarising the key conclusions in this section. 

• There is considerable amount of evidence to suggest that gravitational field equations have the 
same status as, say, the equations of fluid mechanics. They describe the macroscopic, thermody¬ 
namic, limit of an underlying statistical mechanics of the microscopic degrees of freedom of the 
spacetime. The macro and micro descriptions are connected through the heat density Ts of the 
spacetime. Here, the temperature T arises from the interpretation of the null surfaces as local 
Rindler horizons. The entropy density is a phenomenological input, the form of which determines 
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the theory. For a very wide class of theories, it can be defined in terms of a function F(R^,6t) 
built from (2,2) curvature tensor R"[ and the Kronecker deltas, as: 

s = -^ ab ^ p cd=^ V.J3 = o (120) 

• Given a Vector field v a , one can construct three currents: (a) the Noether current J a (v), (b) 
the gravitational momentum P a {v) and (c) the reduced gravitational momentum P a {v). Inter¬ 
estingly enough, one can attribute thermodynamic meaning to these quantities which are usually 
considered to be geometrical. For example, the conserved current </“, associated with the time- 
development vector £“ of the spacetime, leads to a conserved charge (i.e., integral of u a J a (£) 
defined either on a spacelike surface or on a null surface) that can be related to the boundary 
heat density Ts, where T is the Unruh-Davies temperature and s stands for entropy density. 

• One can also define the notion of gravitational momentum P a for all the Lanczos-Lovelock models 
of gravity such that V a (P a + M a ) = 0 (where M a is the momentum density of matter) for all 
observers, leads to the field equation of the Lanczos-Lovelock model. This generalizes a previous 
result for general relativity. 

• The field equations can also be derived from a thermodynamic variational principle, which essen¬ 
tially extremises the total heat density of all the null surfaces in the spacetime. This variational 
principle can be expressed directly in terms of the total gravitational momentum, thereby pro¬ 
viding it with a simple physical interpretation. 

• One can associate with any null surface the two null vector fields £ a ,fc a with l a k a = —1 and 
£ q being the tangent vector to the congruence defining the null surface, as well as the 2-metric 
q a b = gab + (- a kb + £bk a - These structures define three natural projections of the gravitational 
momentum (P a £ a , P a k ai P a q a b), all of which have thermodynamic significance. The first one 
leads to the description of time evolution of the null surface in terms of suitably defined bulk and 
surface degrees of freedom; the second leads to a thermodynamic identity which can be written 
in the form TdS = dE + PdV; the third leads to a Navier-Stokes equation for the transverse 
degrees of freedom on the null surface which can be interpreted as a drift velocity. 

These results again demonstrate that the emergent gravity paradigm enriches our understanding of the 
spacetime dynamics and the structure of null surfaces, by allowing a rich variety of thermodynamic 
backdrops for the geometrical variables. 
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Appendices 

A General Analysis Regarding Null Surfaces 

We will start with a null vector £ a = AW a B, which satisfies the condition £ 2 = 0 only over a single 
surface, which is the null surface under our consideration. Then we obtain 

£ a X/ a £ b = K £ b (121) 


where we have the following expression for k 

k = £ a d a ln A + k ; k=-ht a d a £ 2 (122) 

The last relation defining k can also be written as V a £ 2 = 2 k£ a - The derivation of the result goes 
as follows, let us expand Vb£ 2 in canonical null basis, i.e., V a £ 2 = C£ a + Dk a + Eac Then both 
Ea = e^V a £ 2 and D = —£ a X/ a £ 2 vanishes, since variation of £ 2 along the null surface vanishes. This 
shows that the only non-zero component of V a £ 2 is along £ a . Then it turns out that [43] 

Vi? = Q + K + k (123) 

where, 0 = q ma q m b^ a £ b - Note that the term k enters the picture as £ 2 = 0 only on the null surface. 
With this setup let us now find out R a b£ a £ b in detail, which leads to, 

Rab£ a £ b = £ j (ViVj-r - V,V,/’) 

= Vi {£ j Vjt) - Vj (FVif) - + (V,f ) 2 

= Vi (eVjP - tVjff) - (vT'VA - (Vif ) 2 ) (124) 

However in general, V VjP = k£ 1 is not true, it only holds on the null surface (when £ 2 = 0 everywhere 
this relation is also true everywhere). Since we were doing the calculation for the most general case, 
£ 2 ^ 0 in the above expression we cannot substitute £ : 'V :l £ l = kR, since it appears inside the derivative. 
Thus for the special case when £ 2 = 0 everywhere, we will arrive at the following result 

Rab£ a £ b = -V, (0f) - (vt-'V/ - (V/) 2 ) (125) 

In order to simplify things quiet a bit we will compute the last term (VjPS/jR — (Vjf”) 2 ^ which we 
designate by S. Then we start by calculating the following object on the null surface 

© ah = Q a m qy m r 

= (C + «m + k a £ m ) (S b n + £ b k n + k b £ n ) V m r 
= S7 a £ b + £ a k m \7 m l b + nk a £ b + £ b k n W a £ n + £ a £ b k m k n V rn £ n 
- nk a £ b + k£ a k b - k£ a k b + nk a k b £ 2 

= X/ a £ b + £ a k m V m £ b + l b k n W a £ n + £ a £ b k m k n V m £ n (126) 
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In arriving at the third line we have used the following results: 4V a 4 = k£ 1 and £ a V b £ a = k£ b . Then 
we can reverse the above equation leading to 

vj b = e ab - £ a k m v m 4 - t b k n v a r - 44 k m k n s7 m r ( 127 ) 

From the above equation we can derive two very important identity: 

(v a 4) (v a 4) = e ab e ab + £ a £ b k m v m £ b k n \/ a £ n + £j b k n v a rk m v m i b 
= e ab e ab + 2 ( kk m i m ) ( ni n k n ) 

= Q ab Q ab + 2 kR (128) 

In the same spirit we will arrive at 

(v a 4) (v b r) = e ab e ab + 44 k n v a rk m x7 b e m + 4 ^ m v m rr v „4 

= QabQ ab + ( hk m £m ) (Kk n £n) + 44.A") 444 m ) 

= 0 ab 0 ab + K 2 + K 2 (129) 

The extrinsic curvature for null surfaces, i.e., 0 ab can be given a very natural interpretation. This 
essentially follows from [43]. There the expression for 0 ab in terms of Lie variation of q ab along the null 
generator 4 was obtained as, 

Qab = \q™q%£iq mn (130) 

Now expanding out the Lie derivative term we obtain, 

— £ diQmn T T qandm£ (131) 

Which on being substituted in Eq. (130) immediately leads to, 

Qab = \q™q?£ l diq m n + ^q a iq b d n £ l + ^q b iq^d m £ l (132) 

Now on the null surface q ab = qAB as the only non-zero component. Hence the above equation can be 

written as, 

Qab = 0AB = 2~d\ qAB 2 qAC ® B ^ C + 2 qBC ®A£ C (133) 

On the null surface q%£ b = 0, which in this coordinate system leads to i A = 0 on the null surface. Since 
8a £ 2 represent derivatives on the null surface it also vanishes. If 4 = 0 everywhere, then also t A is 
identically zero everywhere. Hence we have 

Qab = l~k qAB (134) 

There is another way to get this result. If e A are the basis vectors on the null surface and if £ a ,e A 
forms coordinate basis vectors, then qAB = qab^A e B a sca l ar under 4-dimensional coordinate trans¬ 
formation. This immediately leads to the previous expression. For more discussions along identical 
lines see [43]. 
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Now the expression for the quantity S can be obtained as 


5 = V,/'V,(' - (Vif ) 2 

= e ab e ab + k 2 + k 2 - (e + « + it) 2 

= (0 a h© ab - B 2 ) — 20 (k + k) — 2kk (135) 

Using the general expression for R ab £ a £ b we obtain the following form: 

R ab re b = Vi - eVjV) - s 

= Vi (PVjP - UV,-4) - (QabQ ab - e 2 ) + 20 (k + ft) + 2 kk (136) 

For the situation where, l 2 = 0 everywhere we finally arrive at the following simplified expression 
# a6 r£ 6 = - (0 a fc0 ab - 0 2 ) + 20 k + Vi (kU - [0 + k] t) 

= - (0afc0 ab - 0 2 ) + 0K - ~4 (v^0) (137) 

v'q aA 

Let us now try to derive the Raychaudhuri equation starting from the basic properties of null surfaces. 
We start with the following result 

rv Q (v c 40 = rv a v c £ d 

= i?d6ac« 6 +rV c V a £ d 

= V c (fV a 4) - V Q ^V c r - R 6dac £ b r (138) 

Then contraction of the indices c, d leads to the following result 

rv Q (VcO = v c (fv a f) - v a 4v b r - i? a6 « b ( 139 ) 

Otherwise we can rewrite it in a different manner which exactly coincide with Eq. (124). On using 
Eq. (129) and the decomposition: Q ab = (1/2 )0q ab + a ab + u] ab we arrive at 

rv a (V c £ c ) - v c (rv a 4) = -© a6 © ab - « 2 - h 2 - R ab £ a £ b 

= ^ I© 2 _ a °b aab + u ab u, ab k 2 k 2 R ab £ a £ b (140) 

For the situation where, i 2 = 0 the left hand side is just: d (0 + k) /d\ and the first term on the right 
hand side is dn/dX + k(0 + n) the above equation leads to 

^ = K© + k 2 _ v a 4v b r - R ab r£ b 

dX 

= K0 - 0a60“ b - Rab£ a £ b 

= «0 - i© 2 - <J ab Vab + U ab u ab - R ab £ a £ b (141) 

where to arrive at the last line we have used the following decomposition: 0 Q b = (l/2)Qq ab + u ab + 0Jab- 
This is precisely the one obtained in [12] though in a completely different manner. 
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The next object to consider is the quantity £ a J a {£)- This can be obtained by using the identity for 
Noether current leading to, 


j l a J a {l) = V fc [{t a t - fg ab } 

= v 6 1 


V a A 

A 2 


±*(k-s)]-v 6 (/ 2 ^) 
= [( K - «) ^ (« - «) 2 - 


v b 4 


A ' A 2 
= [(« - «) ^ (« - «) 2 - ( K - «) 


(142) 


This can be written in a slightly modified manner as, 

e a r(£) = Vi [(«-s)^] - (At 2 -ft 2 ) 


rVj (k — ft) + (k — ft) (0 + At + ft) — (k 2 — ft 2 ) 

= — (k — ft) + 0 (k — ft) 

dX 


The above expression can be simplified significantly by noting that 0 = d(ln y/q)/d\. 


4 J a W = TA [(K _ s)V g] 


(143) 

which leads to 

(144) 


Again, we have 


D a [(« - ft) T] = ( g ab + tk b + i b k a ) V a [(« - k) 4] 

= V a [(At - ft) r\ + (i a k b + t b k a ) [(« - k) V a 4 + 4V a (At - ft)] 

= v a [(At — ft) £“] — At (At — ft) — ft (At — ft) — — (At — ft) 

dX 

= (At — ft) (0 + At + ft) — (At 2 — ft 2 ) 


(145) 

(146) 


Thus we arrive at 


4 j“(£) = Ha [(/t - ft) £“] + (At - ft) + (At - ft) (At + ft) - (a t 2 - ft 2 ) 

aX 


D a [(At - ft) t] + — (At - ft) 


From the expression of the Noether current we get, 

4 J a (£) = 2 R ab i a t b + 4 g ij £tN% 

The above equation can be used to write, g ab £^N^ b in terms of At and R a b- For that purpose 


(147) 

(148) 
we use 
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Eq. (1) and insert Eq. (143) leading to 

Lg ij £eN% = i a J a {£) - 2 R ab i a t b 

= (k - ft) + 0 (k - S) - 2R ab £ a t b 

CLA 

= ^ (V9«) - 2i? a6 r£ 6 - -±= A [^(« + ft)] (149) 

The above expression when integrated over the null surface with integration measure, yfqdrxdX and 
then being divided by 167 t leads to 

y^r J d 2 xdXy/q^i a g' i:l £eN“ j + Wd (« + «)] | = J d 2 Xy/qK\l - ^ J d 2 xdX^qR ab £ a £ b 

(150) 

Then on using the held equation: _R ab = 87 t [T ab — (l/ 2 )< 7 ab T] = 87 rT ab and then being substituted in 
Eq. (150) we arrive at the following expression 

y J d 2 Xy/q( y ^- s j \i~ J d 2 xd\y/q T ab i a t b = y^- J d 2 xd\y/q^ a g i:i £ t N^ + -^=-^ [\/9(« + «)] 

0 

where the last equality follows from the fact that £ 2 = 0 on the null surface. The above equation 
be written in a more abstract form as 

it J = [1 £ (£) - 3 1 ^ (£)_ 

- J d 2 xdXyfq T ab ti b - yi- J d 2 xdX-^ [y/q (k + k)] 

As an illustration when £ 2 = 0 everywhere, we have R = 0, then Eq. (152) leads to, 

4- J d\d\ vW^Jvg = \ [£ (£)-/, (t) 


(151) 

can 


(152) 


- J d 2 xd\y/q T ab i a t b 


(153) 


We will now try to obtain an expression for the quantity £ a g ,J £(N-‘ :i independently. For that we start 
with the symmetric and anti-symmetric part of the derivative V a i b such that: 

s ab = w a t b + v fe r, j ab = y a t - v b r 

Then we have the following result: S7 a £ b = (1/2 )(S ab + J ab ), which on being 

V b (V a £ b ) - V a (V b £ b ) = R%£ b 


(154) 

substituted in the identity: 

(155) 
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leads to the following identification: 


g ab £,N c ab = - V 6 (S bc - g bc S ) (156) 

Hence we arrive at the following relation: 

l a g bc £^ c = -£aV b [V a 4 + V b r - 2g ab S7d c ] 

= - v b [rv a 4 + 4 v b r - 24 (v c r)] + vj b v b e a + v„4v°4 - 2 (v c 4) 2 

= - v h [rv a 4 + 4v b r - 24 (v c r)] + 2& ab e ab + (k + k) 2 - 2 (0 + « + s) 2 

= -v b [rv a 4 + 4v fc r - 24 (v c 4)] + 2 ( & ab e ab - e 2 ) - 40 (« + «) - (« + k) 2 (157) 

For the case 4 = 0 the term within bracket can be written in a simplified manner such that Lie 
derivative term gets simplified leading to: 

i a g bc £iN bc = 2 ( Q ab Q ab - 0 2 ) - 40k - k 2 + V 6 [(20 + k) 4] 

= 2 (e ab Q ab - 0 2 ) - 40k -k 2 + (20 + k) (0 + k) + 4 (20 + k) 

CLA 

= 2 (0 a b0 Qb - 0 2 ) - 0K + 4 -k + 2 JL (^0) (158) 

v 7 dX y/q dX 

If the null generator is affinely parametrized then k = 0 and Eq. (158) reduces to: 

£ a g bc £ e N bc = 2 (0 o6 0 ab - 0 2 ) + 44 (v 7 ?©) (159) 

' y/qdX 

While for the null generator l a in GNC we have (see Appendix B Eq. (200)): 

= 2 (e„„e*‘ + + (V?K) (wo) 

Through the above analysis we have obtained expressions for t a J a {t), R ab £ a £ b and t a g bc £i.N° ic . 

It turns out from the above analysis that 0 a b0 ab — 0 2 can be given a more physical meaning by 
considering Lie variation of gravitational momentum. This can be obtained by considering variation of 
the gravitational momentum first: 

q ab 6U ab = q ab S [^q(@ ab ~Qq ab )] 

= qabVvSQ ab - 2^0 - v 7 ? QqabSq ab - 

= Vqq ab SQ ab - 2y/q5Q + QS^q (161) 

Now specializing to Lie variation we arrive at: 

-q a b£eR ab = -Q£ty/q - sjqq ab £i® ab + 2y/q£ e G 

= -y/q£ e & + qQ ab £ e q ab - &£e^q + 2y/q£ e G 
= 2^/q ( Q ab Q ab - 0 2 ) + £ t {y/qG) 

= y/ql a g bc £iN bc - (162) 
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where in the last line we have used Eq. (159). Here the quantities 0 ab and 0 can be defined as: O ab = 
( \/2)£tq a b and 0 = i^ln y/q. In GNC parametrization: (1/2 )£(.q ab = (1/2 )d u q a b and £ehiy/q = 
d u In yjq. Thus the Lie variation of gravitational momentum for affine parametrization is directly 
related to T >, i.e. to (0 a f,0 ab — 0 2 )- 

For non-affine parametrization the gravitational momentum associated with null surfaces can be 
taken as n ah = y/q(Q ab — (0 + n)q ab ). Then we readily arrive at the Lie variation expression: 

-q ab £eYl ab = 2 sfq (9 ab @ab ~ 0 2 ) + £( (v 7 ? 0 ) + %V<i £ e K (163) 

Since k is a scalar Lie variation term can be written as: £^k = dn/dX, where A is the parameter along 
i a . If we consider the null generator i a from GNC then we arrive at (see Eq. (160)): 

-q ab £dl ab = ^£ a g bc £eK c - ^ (y/qn) (164) 

while if we have defined the conjugate momenta n ab as n ab = y/q(@ ab — Qq ab ), then the above relation 
could have been written as, 


-q ab £ e Yl ab = ^ a g bc £ t NZ c - ^ + 2k A (^) (165) 

B Derivation of Various Expressions Used in Text 

This appendix will contain derivations of most of the results that we have used in the main text. The 
derivations will be arranged in the same order as that in the main text. First we will present derivations 
related to the Navier-Stokes equation and then we will present the requisite derivations of subsequent 
sections. 

B.l Derivation Regarding Navier-Stokes Equation 

The first thing to compute is the Lie derivative of the object N£ b . This can be obtained starting from 
the first principle, i.e., using expression for N£ b in terms of Y bc and then using Lie variation of the 
connection. This immediately leads to: 

f M a — n ad f r e 4- n ad f r e 

°° v iy bc V be cd ' ^ce bd 

= Qbi (V c v d u e + R e dmc v m ) + Q ad (VfcVdU 6 + R e dmb v m ) 

= i (c 5%V c V d v d + S a c V b V d v d ) - i (V b V c v a + V c V b v a ) 

-\{R a bmc + R a cmb )v m (166) 

In the above expression the second term in the last line can be written as: 

(V 6 V c u a + V c V b u a ) = 2d b d c v a + 2 T a bd d c v d + 2 T a cd d b v d - 2 Y d bc d d v a 

+ v d ( d b T a cd + d c Y a bd ) - 2Y d bc Y a de v e + Y a bd Y d ce v e + Y a cd Y d be v e (167) 
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In order to compute the Lie variation of N£ b along the transverse direction we need the two objects 
£(N A r and £(N BC . For the evaluation of £t.N A r the following identities will be useful: 


ur 

nA , OT^ 


(V fe V c r + V c V b £ a )i = 2 d u d r £ A + 2 V A d d r £ d + 2 Y A d d u £ d - 2 T d ur d d £ A 

+ v d (d u v A d + d r r A d ) - 2Y d ur v A e t + r A d r d re v e + v A d v d ue t 
= 2d u p A + 4aT A + 2/3 b T a b - 2p A T r ur + d u T A u + d r T 


A 

uu 


_ opd pA _i_ pA pd _i_ pA pd 

urdu ' ^ ud*~ ur ' rd*~ uu 


(168) 


and 


(v fe v c r + v c v 5 rp 


+ \{R\ mc + R a cmb ) Hi = 2d u ^A + 4ar A r + 20 B r A B - 20 A K 

o -pA _|_ o pA _ 9F'^ pA _i_ P^ P^ _i_ P^ P^ _i_ 

l U U L ru ~r U r L uu ur L 1 U( J^L ur ~T 1 r( ^L uu “p -it 

= 25 u /3 a + 2(3 B r i B + 25 u rt 
= d u /3 A + p B q AC d u q B c 


(169) 


While for £(N BC we have: 


(v h v c r + v c v b r) A c + [(R a bmc + R a cmb ) H BC = - 2 /Ar bc + 9 B r A c + d c r A B - 2 T d BC r A d 


+ ^B<AuC + ^Cd^t,B + BC — dc^Bu + 

pA pd , o fA o pA I pA pd pA pd 
— 1 CcP itB °u*- BC ~ Cu “I" 1 nd 1 BC ~ 1 Bd l uC 

= —2f5 A T r BC + 2d u f BC 


(170) 


These are the expressions used to get expressions in Section 7.1. From the vector Yl a given in Eq. (96) 
we can calculate the Lie variation along £ a leading to, 


and equivalently, 


£ t n n = £ m d m n n + n m d n r 
= fo, ^Pa/ 3 a , \d u p A 


£tSi n = £ m d m n n - n m d m r 


= (o,o, -d u p 


(171) 


(172) 


Also, 


D m O™ = d B Q B A + d c hi ^& C A - Y C AB Q B 

= \d B ( q BC d u q A c ) + ^ q CD d c In ^fqd u q AD - \q BD d u q CD £ C A B (173) 

Using these results we finally obtain, 

q n a £(Si n + D m ep + en a - o a (e + a ) = l -d u e A + d B Q q BC d u q A c ) 

+ i q CD d u qADdc In v 7 ? - + du In \fq\pA - d A d u In y/q- d A a (174) 


37 


In raising the free index of the above equation the following identities can be useful 

- q CD q AB d u q BC d D In y/q = d u (q AD d D In y/q) - q AD d u d D In y/q (175a) 

d u ( q AD d D In y/q) - q AB d u q CF d c qBF ~ q AB d D ( q CD d u q B c ) 

= ~d u (q BC £+ q CF d u q AB d c qBF - q AB d D q CD d u q BC 
= -d u ( q BC f A c ) (175b) 

q AB q CD d u q E D^ B c = f p C d u q CF - q AB d c q B Fd u q CF (175c) 

Moreover we also have: 

Rabt a q b c = RuA = G u a = ^d u /3 A - d A a + ^q BC d u d B q CA + ^ d B q BC d u q AC - d u d A hi y/q 

+ ^ P A d u In y/q + i q BC d u q AC d B In y/q - ]^q BD d u q C D^ C AB (176) 

It can be checked that this expression coincides exactly with Eq. (104) in Section 7.1 as it should. 

To bring out the physics associated with Noether current and its various projections we compute 
the Noether potential and hence the Noether current completely in GNC for the vector To start 
with we provide all the components of the tensor V a £b, which are: 

/Ja(,b) uu = ~rd u a; (V a £ b ) ur = a + rd r a; (V a 6) ru = -a - rd r a (177) 

(Va£,b) uA = rd A a-rdul3 A ; (V a &) Au = -rd A a\ (Va€b) rA = ~\&A ~ \rd r P A (178) 

(^aCb)rr = a£b) Ar = 7/pA + — rd r /3 A ', 

(V a(,b) AB = l;d u q AB + ir ( d A /3 B - d B fi A ); (V Q &) Bj4 = ^ d u q AB - ^r (d A /3 B - d B /3 A ) (179) 


Then components of Noether potential J a b = V a £b — V&£ a have the following expression: 

Juu = 0; Jur =2a + 2rd r a; J uA = 2 rd A a - rd u /3 A (180) 

J r A = -/3a - rd r /3 A ; J AB = r ( d A /3 B - d B j3 A ) (181) 

The upper components of Noether potential can be obtained as 

juu = 0 . jur = _ 2a _ 2r g rC( _ r p A pA _ r 2pAg r p A ( 18 2) 

JUA = _pA _ rq AB dr p B . jrr = ^ pA p B _ Q bM ( 183 ) 

J rA = 2rq AB d B a - rq AB d u p B - 2r 2 aq AB d r p B - r 3 p 2 q AB d r p B 

- r 2 q AB p c (d B Pc ~ d c pB) + 2r 2 p A d r a - r 3 / 3 A p B d r p B (184) 

J AB = - r p A q BC (/3c + rdrPc) + rp B q AC (Pc + rd r p c ) + rq AC q BD (d c p D ~ d D p c ) (185) 


Using the above components of Noether potential the components of Noether current can be obtained 


38 


as 


(186) 

(187) 

(188) 


J u (0 = -4 d r a - p 2 - 2ad r In Jq - -L d A (^ A ) 

J r (£) = 2a5 M In ,/g + 29 u a 
J A (^) = ( v / 9/3 A ) + q AB d u p B - 2q AB d B a 

Note that k a J a {£,) = — J“(^), g b J b (£) = J A (£) and finally l a J a (C) = J r (£)- As we will see & h °f them 
matches with our desired expressions. Also in the stationary limit we have d u a = 8 u Pa = d u qAB = 0, 
which in particular tells us that J r = 0. Hence in the static limit Noether current is on the null surface 
since its component along k a (which is —£ a J a (£,)) vanishes. Also in this case we have: 


(VfcVcf +v c v 6 nt + [(# 

= 2 a,r 


bmc 

A 


, rja \ tmi-4 rs pA , n y/ 
' cmb) s J ur u u L ru i u r 1 u 

= -d u p A 


■pa "pAi 

ur*- du 


pA -pd 
*- rd*- uu 


+ R / 


(189) 


as well as, 

(V b V c r + V c V 6 C°)bc + K Ra bmc + R a cmb ) flflc = d B^uC + 9 C T A B - 2 T d BC T A d 

+ ^Bd^tc + ^CcAIiB + — dc^ Bu + ^ud T BC 

t^A j^d i cj f^A q t-'A i pA p'd t-'A -pd 
~ 1 Cd 1 uB “T Out BC ~ °B L Cn 1 ud 1 BC ~ 1 Bd L uC 

= 2d u t BC 

Using these two results we arrive at: 


£pN a = id, 


£sN A c = -6 A d c e 


1 


6£d B e - d u v A c 


These are the expressions used in Section 7.1. 


(190) 


(191) 

(192) 


B.2 Derivation Regarding Spacetime Evolution 

We need to consider the object lag 1 -’ £(.Nij in GNC. This in turn requires us to obtain expressions 
for £eN£ r and £g_N T AB . Then using the identity for Lie variation of N° b we can obtain both the Lie 
variations. For that purpose we have: 


1 

2 


(s%v c v d i d + 5“VbVdi d Y = \d u e + d u 

\ / ur Z 


= 2.d u a 


- - (R 


= 0 


fv b v c r + v c v b r) ? 

0 

fv 6 V c r + V c V b r) = ad u q AB - \q CD d u q A cd u qBD 

\ ) AB 2 

K R “ 


I c>a 
bmc ' cmb 


r 1 ]_ \ r -< 

= --Oid u qAB + -zduQAB ~ ~;Q. CD d u qAc9uQBD 
ab 2 2 4 


(193) 

(194) 

(195) 

(196) 

(197) 
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This immediately leads to 


£tKr = \%^y/Q (198) 

£(.Nab = ~ a d u qAB + 2^^ ab (199) 

Combining all the pieces and using the results 0 = d u In yfq and Qab = 0-/2)d u qAB, which is the only 
non-zero component of Q a b, [43] we finally obtain 

£ a g ij £eNf. = 2 £ t N r ur + q AB £ e N r AB 

= -2ctd u In y/q + 2d 2 u In y/q - ^ d u qABd u q AB 

= 23.0 + 2 (e„»e- - e a ) + //// - (Vi*) (200) 

which can also be obtained from a completely different viewpoint. For sake of completeness we will 
illustrate the alternative methods as well. For the null vector £ a in the adapted GNC system we have: 

(rV c r)“ =a + r/3 2 + 0(r 2 )- (£ c V c £ a ) r = rd u a + 2 ra 2 + 0(r 2 ) 

(t c S7 c t a ) A = ra/3 A + rq CA d c a + 0(r 2 ) (201) 

Hence on r = 0 surface, we have k = a, as well as, R, = —(1/2 )k a 'S7 a £ 2 = a. Now we will use the 
Raychaudhuri equation to get R a b£ a £ b and hence the Lie variation term. In this case we have, du = dX , 
thus Raychaudhuri equation reduces to the following form (see Eq. (139)) 

rVa (0 + 2a) = V c (fv.f) - V a 4V b r - Rabtt (202) 

Where, the 0 + 2a term comes from Then we have, 

v c (rv a £ c ) = d c (rvj c ) + rv a £ c s c in ^ 

= 2 a 2 + ad u In yfq + 2d u a (203) 

Thus non zero components of B a b = V a £b are as follows: 

B ur — Cfc; B rA — 2 @Ai B A C = ~^^uQ_AC (204) 

From which it can be easily derived that, B a bB ba = 2a 2 — (1 /£)d u qABd u q AB ■ Thus we obtain 

Rab£ a £ b = -d u Q + 2a 2 + Qa - B ab B ba 

= a© - \q AB d 2 u q AB - jd u q A Bd u q AB 

= a©- l -d 2 u ^/q + (d u In y/q) 2 - Q ab Q ab (205) 

where <d a b has the only non-zero component to be, Qab = (1/ 2)d u q A B ■ For the GNC null normal £ a , 
the Noether current vanishes, such that Lie variation of N£ c turns out to have the following expression 

= -2 R ab ri b 

= 2 d u a + 2 (QabQ ab - © 2 ) + -^ ^ (y/qa) (206) 
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The components of S a b = V„,4 + V&4 In GNC are as follows: 


S uu = 2 rd u a - 4ra 2 + 0(r 2 )\ S ur = 2a + 2 rd r a + r/3 2 + 0(r 2 ) 

S U A = -rP B d u qAB + 2rdAOi - 2 raf^A + 0(r 2 )- S rr = 0 
SrA = Pa + rd r /3 A ~ r(3 c d r q c a + 0(r 2 ) 

Sab = d u q A B + 2rad r q A B + r ( D A /3 B + D b /3a) + 0(r 2 ) (207) 

Thus the trace at r = 0 leads to: 5 = 4a + 2d u In yfq. Thus we arrive at the following expression (see 
Eq. (156) of Appendix A) 

Lg ij £iN“ = 2d u (0 + 2a) - d b S rb - T£ c S bc - S rc 8 c In yfq (208) 

Then the upper components of S ab necessary for the above computation are the followings: 

S ur = S ur + r/3 A S rA = 2a + 2 rd r a + 2r/3 2 + 0(r 2 ) 

S rr = 2rd u a + 4ra 2 + 0(r 2 ) 

S rA = 4ar/3 A + 2rq AB dACt — 2ra/3 A + 0(r 2 ) (209) 

The mixed components leads to nothing new so we have not presented them. From Eq. (208) the 
expression for Lie derivative can be obtained as: 

4 g ij £eN% = 2d u (0 + 2a) - d u S ru - d r S rr - d A S rA + 4a 2 + 20 ab 0 ab - 2a0 

= 2 (0 ab 0 ab - 0 2 ) + -p-rr (v 7 ?©) _ 2a0 + 4 d u a + 4a 2 - 4<9 u a - 4a 2 

y/q d\ 

= 2 (0 a6 0 ab - 0 2 ) + 4= 4 t (V?®) - 2 «© ( 21 °) 

y/q aA 

which exactly matches with Eq. (206). The same can be ascertained for Eq. (205) by computing R ab £ a £ b 
on the null surface i.e. in the r —> 0 limit, directly leading to: 

d na nb _ td _ ft pa _ o pa . pa p6 _ pa p6 

Z- £^uu G'a 1 uu ua ' 1 uu 1 - ab ubua 

= duKu + dr Ku + dAT A u - dl In yfq + T u uu d u In yfq- T a ub T b ua 

= 2a 2 - d 2 u ln Vg + ad n in fq ~ Kb^iu ~ Kb^ur - £ub rb uA 

= -dlhiy/q + aduhiy/q- 0 ab 0 ab (211) 

which under some manipulations will match exactly with Eq. (205). Then in GNC we obtain in identical 
fashion, the following expression for heat density, 

5 = vrfjVW - (v.,r) 2 

= (2a 2 - (1/4 )d u q A Bd u q AB ) - (0 + 2a) 2 

= -2a 2 - 4a0 - 0 2 + 0 ah 0 ab (212) 

This on integration over the null surface leads to, 

-5- J dud 2 x yfq S = J dud 2 Xy/q(Qab® ab — 0 2 ) — J dudfxyfqa 2 — 4 J d 2 xTds (213) 
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Let us now write the integral form of R a b£ a l b , for that we note the integration measure to be dud 2 Xy/q. 
Thus on integration with proper measure and (1/87 t) factor leads to 


h. 


1 dA i 


dud z Xy/qR a bt a t = / dud 2 Xy/qD - ^ ^ 


+ / d 2 xTs\ 2 — / d 2 xsdT 


(214) 


which can be written in a slightly modified manner as: 


1 


1 


— / dud 2 xAqR a b £ a £“ = —— / dud 2 Xy/qD — 

' o7T , ' 


1 dM L 


87T 


87t dA 


d 2 xTds 


(215) 


Also the Lie variation term (with all the surface contributions kept) on being integrated over the null 
surface we obtain 


1 

167T 


dud 2 Xy/q x i a g^ £(.N^ 


1 

87T 


dud 2 Xy/qD + 



X Tds + — 
87T 


dud 2 Xy/qD 


1 dA± 2 _ 
87t dA l 
1 dA± 2 

+ 8n d\ l 



(216) 


To calculate Lie variation for we need to calculate V a £b + Vb£ a = S a b- This tensor has the following 
components: 


Ruu — 2r5 u G:, Rur — 0 

SuA — Rrr — 0 

RrA = o. Rab = d u QAB (217) 


Thus in the null limit obtained from the relation: r —> 0 we arrive at the result that all the components 
of R a b vanishes except for the Rab components. If we want to satisfy the Killing condition for £“ on 
the null surface we would require d u qAB = 0. From the above relations it is clear that V a £ a = 0. 
Moreover we also have, 


K = -hCVae = -I= Ku = a (218a) 

k = ~k b V b £ 2 = \d r {-2ra) = ~a (218b) 


which shows that for f a , k = k. Thus even without the condition d u qAB = 0, we arrive at the relation 
k = —k = a. Moreover Lie variation of N£ c along £ a can be obtained by computing the following 
objects: 


1(«V C W{'‘ 

(v,v c r 

— — ( ]V l 

2 bmc 

(v>Vrf“ 

[_ if pa 

2 y 1X bmc 


VcVbCY =-2d u a 

/ ur 

b)r] r =o 

\T ^ 

: V b£, a ) = —Oid u qAB ~ -q CD d u qAcd u qBD 

/ AB 2 

-n CD r 


f? a 

cmb 

-v, 


R 1 


\ I r \ \ \ 

l cmb)C l \ AB = - g ad udAB + 2 d udAB - -gt 


(219) 

( 220 ) 
( 221 ) 
( 222 ) 
(223) 
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which can be used to obtain the Lie variation term associated with £ a as, 


e a g lJ £zN“ = 2d u a + 2 (0 a6 0 ab - 0 2 ) + — dl^ 

V9 

= -^pd u (ay/q) + tag 13 £(N^ 

y/q 


(224) 


Then using the momentum n ab = y/q[Q ab — q ab (Q + k)] conjugate to the induced metric q a b from 
Eq. (163) we immediately arrive at, 


a 


-q ab £ i n a0 = y/jl a g l3 £ i N“- dy 
These expressions are used to obtain Eq. (86). Also the variational principles in this context are: 


(225) 


Qi = j d\d 2 xyfi (-±-R ab i a i b + T ab re b ^j 


= J d\d 2 Xy/^[^-V + T ab e a £ b 


1 


— / d 2 x Tds + 


1 dA l 


Q 2 = J d\d 2 Xy/q 

= J d\d 2 Xy/q 
These are the expressions used in Section 7.3. 


-^-l a g' 3 £tN? j+ T ab n b 


2 ~ ^\^-V + T ab £ a £ b 
L 87 r 


+ / d 2 x sdT + 


8n dX 


1 dA± 


8ir dX 


(226a) 


(226b) 


References 

[1] J. D. Bekenstein, “Black holes and entropy,” Phys.Rev. D7 (1973) 2333-2346. 

[2] J. D. Bekenstein, “Generalized second law of thermodynamics in black hole physics,” 

Phys.Rev. D9 (1974) 3292-3300. 

[3] S. Hawking, “Particle Creation by Black Holes,” Commun.Math.Phys. 43 (1975) 199-220. 

[4] S. Hawking, “Black Holes and Thermodynamics,” Phys.Rev. D13 (1976) 191 T97. 

[5] P. Davies, “Scalar particle production in Schwarzschild and Rindler metrics,” 

J.Phys. A8 (1975) 609-616. 

[6] W. Unruh, “Notes on black hole evaporation,” Phys.Rev. D14 (1976) 870. 

[7] G. W. Gibbons and S. W. Hawking, “Cosmological Event Horizons, Thermodynamics, and 
Particle Creation,” Phys. Rev. D15 (1977) 2738-2751. 

[8] V. Iyer and R. M. Wald, “Some properties of Noether charge and a proposal for dynamical black 
hole entropy,” Phys.Rev. D50 (1994) 846-864, arXiv: gr-qc/9403028 [gr-qc] . 


43 













[9] T. Padmanabhan, “Thermodynamical Aspects of Gravity: New insights,” 

Rept. Prog. Phys. 73 (2010) 046901, arXiv:0911.5004 [gr-qc]. 

[10] T. Padmanabhan, “Gravity and the thermodynamics of horizons,” 

Phys.Rept. 406 (2005) 49-125, arXiv:gr-qc/0311036 [gr-qc]. 

[11] T. Padmanabhan, “Emergent Gravity Paradigm: Recent Progress,” 

Mod. Phys. Lett. A30 no. 03n04, (2015) 1540007, arXiv: 1410.6285 [gr-qc]. 

[12] T.Padmanabhan, Gravitation: Foundations and Frontiers. Cambridge University Press, 
Cambridge, UK, 2010. 

[13] T. Padmanabhan and D. Kothawala, “Lanczos-Lovelock models of gravity,” 

Phys.Rept. 531 (2013) 115 171, arXiv: 1302.2151 [gr-qc]. 

[14] T. Padmanabhan and A. Paranjape, “Entropy of null surfaces and dynamics of spacetime,” 

Phys.Rev. D75 (2007) 064004, arXiv:gr-qc/0701003 [gr-qc], 

[15] T. Padmanabhan, “Classical and quantum thermodynamics of horizons in spherically symmetric 
space-times,” Class. Quant. Grav. 19 (2002) 5387-5408, arXiv: gr-qc/0204019 [gr-qc], 

[16] T. Padmanabhan, “Dark energy and gravity,” Gen.Rel. Grav. 40 (2008) 529-564, 
arXiv:0705.2533 [gr-qc]. 

[17] T. Padmanabhan, “Equipartition of energy in the horizon degrees of freedom and the emergence 
of gravity,” Mod.Phys.Lett. A25 (2010) 1129-1136, arXiv:0912.3165 [gr-qc]. 

[18] T. Padmanabhan, “Surface Density of Spacetime Degrees of Freedom from Equipartition Law in 
theories of Gravity,” Phys.Rev. D81 (2010) 124040, arXiv: 1003.5665 [gr-qc], 

[19] T. Padmanabhan, “General Relativity from a Thermodynamic Perspective,” Gen.Rel. Grav. 46 
(2014) 1673, arXiv:1312.3253 [gr-qc], 

[20] S. Chakraborty and T. Padmanabhan, “Evolution of Spacetime arises due to the departure from 
Holographic Equipartition in all Lanczos-Lovclock Theories of Gravity,” 

Phys.Rev. D90 no. 12, (2014) 124017, arXiv: 1408.4679 [gr-qc]. 

[21] K. Parattu, B. R. Majhi, and T. Padmanabhan, “Structure of the gravitational action and its 
relation with horizon thermodynamics and emergent gravity paradigm,” 

Phys. Rev. D 87 (Jun, 2013) 124011, arXiv:gr-qc/1303.1535 [gr-qc], 
http://link.aps.org/doi/10.1103/PhysRevD.87.124011. 

[22] S. Chakraborty and T. Padmanabhan, “Geometrical variables with direct thermodynamic 
significance in Lanczos-Lovelock gravity,” Phys.Rev. D90 no. 8, (2014) 084021, 

arXiv:1408.4791 [gr-qc]. 

[23] T. Padmanabhan, “Lessons from Classical Gravity about the Quantum Structure of Spacetime,” 
arXiv:1012.4476 [gr-qc]. 

[24] T. Padmanabhan, “Entropy density of spacetime and the Navier-Stokes fluid dynamics of null 
surfaces,” Phys.Rev. D83 (2011) 044048, arXiv: 1012.0119 [gr-qc]. 


44 


[25] S. Kolekar and T. Padmanabhan, “Action principle for the Fluid-Gravity correspondence and 
emergent gravity,” Phys.Rev. D85 (2012) 024004, arXiv: 1109.5353 [gr-qc]. 

[26] T. Damour, “Quclques proprietes mecaniques, electromagnetiques, thermodynamiques et 
quantiques des trous noirs (available at http://www.ihes.fr/damour/Articles/),” Th‘ese de 
doctorat dEtat, Universite Paris 6 . (1979) . http://www.ihes.fr/damour/Articles/. 

[27] T. Damour, “Surface effects in black hole physics,” Proceedings of the Second Marcel Grossmann 
Meeting on General Relativity (1982) . 

[28] R.-G. Cai and S. P. Kim, “First law of thermodynamics and Friedmann equations of 
Friedmann-Robertson-Walker universe,” JHEP 0502 (2005) 050, 

arXiv:hep-th/0501055 [hep-th]. 

[29] M. Akbar and R.-G. Cai, “Friedmann equations of FRW universe in scalar-tensor gravity, f(R) 
gravity and first law of thermodynamics,” Phys.Lett. B635 (2006) 7-10, 

arXiv:hep-th/0602156 [hep-th]. 

[30] D. Kothawala, S. Sarkar, and T. Padmanabhan, “Einstein’s equations as a thermodynamic 
identity: The Cases of stationary axisymmetric horizons and evolving spherically symmetric 
horizons,” Phys.Lett. B652 (2007) 338-342, arXiv:gr-qc/0701002 [gr-qc], 

[31] A. Paranjape, S. Sarkar, and T. Padmanabhan, “Thermodynamic route to field equations in 
Lancos-Lovelock gravity,” Phys.Rev. D74 (2006) 104015, arXiv:hep-th/0607240 [hep-th], 

[32] D. Kothawala and T. Padmanabhan, “Thermodynamic structure of Lanczos-Lovelock field 
equations from near-horizon symmetries,” Phys. Rev. D79 (2009) 104020, 

arXiv:0904.0215 [gr-qc]. 

[33] S. Chakraborty, K. Parattu, and T. Padmanabhan, “Gravitational Field equations near an 
Arbitrary Null Surface expressed as a Thermodynamic Identity,” arXiv : 1505.05297 [gr-qc]. 

[34] S. Chakraborty, “Lanczos-Lovelock gravity from a thermodynamic perspective,” 
arXiv:1505.07272 [gr-qc]. 

[35] T. Jacobson, “Thermodynamics of space-time: The Einstein equation of state,” 

Phys.Rev.Lett. 75 (1995) 1260-1263, arXiv : gr-qc/9504004 [gr-qc], 

[36] T. Padmanabhan, “Momentum density of spacetime and the gravitational dynamics,” 
arXiv:1506.03814 [gr-qc]. 

[37] D. Lovelock, “The Einstein tensor and its generalizations,” J. Math. Phys. 12 (1971) 498-501. 

[38] C. Lanczos, “Electricity as a natural property of Riemannian geometry,” 

Rev. Mod. Phys. 39 (1932) 716-736. 

[39] C. Lanczos, “A Remarkable property of the Riemann-Christoffel tensor in four dimensions,” 
Annals Math. 39 (1938) 842-850. 

[40] T. Padmanabhan, “Some aspects of field equations in generalised theories of gravity,” 

Phys.Rev. D84 (2011) 124041, arXiv: 1109.3846 [gr-qc]. 


45 


[41] V. Moncrief and J. Isenberg, “Symmetries of cosmological cauchy horizons,” 

Communications in Mathematical Physics 89 no. 3, (1983) 387-413. 

http://dx.doi.org/10.1007/BF01214662. 

[42] E. M. Morales, “On a Second Law of Black Hole Mechanics in a Higher Derivative Theory of 
Gravity,” available at 

http://www.theorie.physik.uni-goettingen.de/forschung/qft/theses/dipl/Morfa-Morales.pdf (2008) 

[43] K. Parattu, S. Chakraborty, B. R. Majhi, and T. Padmanabhan, “Null Surfaces: Counter-term 
for the Action Principle and the Characterization of the Gravitational Degrees of Freedom,” 
arXiv:1501.01053 [gr-qc] . 

[44] S. A. Hayward, “Unified first law of black hole dynamics and relativistic thermodynamics,” 
Class.Quant.Grav. 15 (1998) 3147-3162, arXiv:gr-qc/9710089 [gr-qc]. 

[45] D. Kothawala, “The thermodynamic structure of Einstein tensor,” 

Phys.Rev. D83 (2011) 024026, arXiv: 1010.2207 [gr-qc]. 


46 


